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On Transient Solutions of the “Baffled Piston” 
Problem 


F. Oberhettinger ' 


luced by the movement of a piston n infinite rigid wall for 


1. Introduction 


The euse of the time harmon movement of a piston membrane im an mfinite rigid wall 
“baffle’’) can easily be generalized to the case when the motion of the piston Is not periodic 
but an arbitrary function of time. Such transient solutions have become of considerable in- 
terest in recent times (for a detailed treatment of the propagation of such sound pulses, see 5]? . 
The prom edure for the case treated here is the same as used elsewhere [7], 1.e., the Green’s funce- 
tion for the exponential decay case (modified wave equation Al vl 0,4 ik) is used to 
obtain the solution for the pulse problem The aceustic field velocity potential for the time 
harmonic movement of the piston include representations given by Bouw kamp 1], King [6], 
and Wells and Leitner [9] The first of these contributions gives the solution in the form of a 
series expression while the second and third involve integral representations that are obtained 

75] 


The se repre sentations can be used to treat the general case of an arbitrarv move- 


using integral transform methods Hankel transform |4 p 73] and Lebedev transform +, |] 


respectively 
ment of the piston In view of the method to be emploved here, such representations should be 
used for which the inverse Laplace transform of the velocity potential with respett to the purely 
lagimnary Wave parameter ¥+ k can be given Such an expression can be obtained in a direct 


way by regarding each pom of the moving disk us an accustic point source and integrating 


over all pomts of the disk 


2. Exponential Decay and Time Harmonic Solution 


The piston is represented by an infinitely thin cireular disk of radius a located in the z,y- 
plane with its centet at the origin of a three-dimensional Cartesian svstem of coordinates The 
remaining part of the a) plane consists cl a rigid wall The movement of the disk is time 


harmonic of the form 


(1) 


where is the velocity of the disk 
Each point of the disk can be considered as a point source with the velocity potential 


i 
? ” f po’ de'dp’, PQ, 2) 


ke, cis the velocity of sound which represents Green’s function of the wave equation in 
. | | 


free space (see fir. | 
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Upon substituting 


one obtains, after isolation of the time-dependent factor ¢7* for the free-space Green’s func 


tion of 
Au—vy*‘u=—0, t 
the expression 
ve 
f , ; , ~- 
u U] p de dp » » 
477 
. [n2 ,2 9) fan " “ 2 . 
/ ip p zpp COs (l¢ ¢¥ 2 (th 


The transition from & to y by (3) amounts to the transition from a wave problem to an expo 
nential decay problem. At first the exponential decay problem for positive real y is solved and 
the solution for the wave problem can be obtained by returning from y to 4 This, however 
is not necessary here, since the exponential decay solution will be used in order to solve the 
problem of an arbitrary moving disk. If the expression in (5) is integrated over the whole disk 


one obtains observing (6 


, n , se ty ; 
c = |p? +-p’*—2pp' cos (y—¢ ) 
4r ' 


~~ Of 


. exp p- a’? 2 pp COS (eo (ta) ~ p dp de 7 


Here is the ‘‘densitv distribution” over the disk. Since the disk moves as a W hole, this prop 
n : pro} 


erty is a constant; it is determined by the fact that according to (1 
ol 1, for p<a, 
(—) . 
O02 0, for p>a 


This function Ul’ as defined in (7) satisfies 


(a) AU’ y'l 0, 
ol’ : ol ; 

(b) —1 for p<_d, O for p >a at 2<=—0. 
O: 0 

(c) lim 17-0 when p or 


2 





In order to perform the integration in (7), one replaces the integrand in (7) using the well- 


known formula [3, vol. 1, p. 191, 9 
mS A(t 5 ( di é FA >) 


or, upon subs tituting 


\ r J (Ar) (7 pt lai at (9) 


f 


If this formula is used to replace the integrand in (7), one obtains upon interchanging the 


order of integration 


n : "a "as ; 9 , , ae, 
1 ' r(r . ( . { J (rt p 0 2pp cos (o—-e p dp’ de har 


(10) 


To evaluate the Integral in 10 one uses the addition theorem for the Bessel functions [4, 


p 10] 


J rip p Zpp cos l¢ ¥ a eS Tp J Tp COS | Tile eo 


Integrating with respect to ¢’ and observing that onlv the term m=O gives a contribution 


| . . 
7 T(T " ( : Tp { p’ J rp’ )dp’ dr. 


’ 


\ ields 


one obtains 
| . 
/ an | ret. ents Ji(tp) Jy (ra)dr (11) 


Now necordaing Lo l 


p ad: )} p ad 
O \ 0 
therefore by [3, vol. 2, p. 14 
(<= ) | ad re) ad 
“an 
O 2 () p >a 
Hence 
- > 
Che final re presentation becomes therefore 
U=e 2), Jy (rp) J, (ra) dt (12) 
> {) 
Clearly U7 satisfies the conditions a, b, ¢, stated following (8). The restriction -O mn 12 


means ho loss of generality , Since 1t is sufficient to consider the upper half space lf one wishes 


3 











to return to the wave problem one has to replace y by 7k according to (3) and one arrives in 
this case at King’s solution [6]. 
Approximate evaluations of King’s and related integrals for various cases have been 


carried out [2, 6, 8]. 
3. Transient Solutions 


The expression (12) is now used to investigate the corresponding pulse problem. Let the 
0; let the pulse function be 


movement of the piston be represented by a pulse starting at ¢ 
auad- 


g(t), and g(t)=0 for <0. Under the assumption that the pulse function g(f) is such that it 


mits its representation by Laplace’s integral formula 


g(t) (277) ' | || qr ie rdr |e dy, » 


one can represent the expression for the acoustic field generated by the pulse in the form 


(t)=(2mi)~! | | | g(r)e var [Ue 14 


where lL’ is the representation (12) for the ‘exponential decay” field. (The velocity of sound 
has been chosen to be unity.) 
If the order of integration can be interchanged one obtains 


e()= | a(x) | (ni) i Ue ray ar 5 


The integral in the bracket is the inverse Laplace transform of the “exponential decay”’ solution 


with respect to the wave parameter y. If this is known, then the velocity potential P(t) of the 


acoustic field for the case of an arbitrary movement of the piston can be obtained by a further 


integration involving the pulse function g(t) according to (15). The inverse Laplace trans 


form of the integrand in (i2) with regard to y is known [3, vol. 1, p. 248] 


%*e+i1@ 


(29) 1 (r?+-7”) 2 ¢ 2(r2+y? 2 ely 


¢€ t@ 


Jdr(t— 2 ial 
Lo 
according as t>z or t< 2, respectively. 


The acoustic velocity potential due to a pulse function g(t) (“Dirac”’ pulse movement of 


the piston at ‘=0) is by (14) 


©, (t) = 2ri) U,"dy. 17 
Hence, (15) can be written as 
v= | g(r)® D(t—r) dr. IS 
0 
Therefore, by (12) and (16) 
0, t 
® p(t) = an 1 19 
a| J(ra) Jy(rp) Jo[r(t?—27)* |dr, t 
J0 
The integral is known. One has [3, vol. 2, p. 21]: 
(0), 0< y¥< |b—« 
= }? ? ws 1 
Jy (br) Jy (er) J, (xy)dxr (ry) ~' are cos ( conde ), b—e|\<y<b+e 2) 
Jo ; 2be , 
yo! y>b+« 


In order to evaluate (19 identify y with a, 6 with (2 2)'2 and ¢ with p Two different cases 


have to be distinguished 
Case A a> p, 


Case B a 


1.e., the two possible different cases where the projection S ot the poimt of observation P into 
the LY plane lies either inside or outside the disk are treated separately hig 
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The largest and the smallest distances of a point on the circumference of the disk from the 
pom|l of observation P are denoted by Rand R’ respectively. 


R o*-+- (a-+ p)*]*'*, R’ 2+. (qa—p)?|'/?, (21 


One obtains then for ®)(¢ by (19) under consideration of (20) and (21 
kor the CUS A, a p 


p / 0). t 
l, t R’ 
99 
/ ) a- 9 > 
T arceos | > rE : } R t hi 
«~~ 
() f R 
/ fhe CAN¢ B a fs) 
?, (), t R’ 
/* oe oo a P 
1 arceos | D1 : } R t<R 23) 
=p | e 
) f R 








The physical interpretation of these results is obvious The parts in (22 and (23) that 
are different from zero represent the transient acoustic field at a point P due to a “Dirae’ 
pulse movement of the piston for a time which lies between the smallest and the largest distance 
of the point of observation ? from the disk. Outside of this the field is zero Note that the 
velocity of sound was assumed to be unity.) The “Dirac’’ pulse solutions (22) and (23) ean 
now be used to construct the acoustic field in case of an arbitrary time di pendency g(t) of the 


movement of the piston according to (18). One obtains immediately 


For the case A, a™>p: 


"t-R . 
(7 T po ad . 
®(t)—r ‘| d(r) arecos dr al(r)dr 24 
Ji-R ~p|(t—r . 


For the case B, a<p 
f ‘4 / T 0 a . ~ 
q(T arecos 1 aT 2) 
) ; . 
° K - 


The Cuse p 0) (point of observation on the axis of the disk) reduces to 


(tf T 





It has to be remembered that q(t 0 for <0 Therefore no possibly negative part of thi 
interval of integration in (24) to (26) gives a contribution Furthermore, the time parameter / 


in the above formulas has to be replaced by ct. where ¢ is the velocity of sound 
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Special Types of Partitioned Matrices 


Emilie V. Haynsworth” 


ior direct prodcucts is given mn 


1. Introduction The general definitior 


sector 


| rrecent papel L, ‘) (Cloddard and I] scl eet lt Was Show! by thr iuthotr l ©} that riven a 
} discussed the relation ship betwee! fhatrices | Hatin 1 satisiving >) for some tnatrices B and YY 
ind B. of orders nad 7 respect vely which satisfy there @NXISts a tmatrix P which Cal be constructed 


from YX, such that 
1Y—NB 


ior so i } itrin \ of rank () iP? ( ), } 
{) f 


The res ilts ol Gioddard and Sechnemer were vel 
eralized 3} i paper by the authol >} for partitioned 
trices “1 and 6b which corresponding sub- | so that prof the roots of A are the roots of B. Thus, 
matrices (or blocks | nnd i Satisty i! alogo is il or ol B is reduced } rives a red iction for 
relationship to that ven by (1 mula for A Naturally a itrix satisfying l 
' Sutisaes >) with p | 
1,,.\ Ab , - | theorem 2 of this paper a similar reduction 
. lor ila Lo those for matrices satisiving | and o 
vhere 1), 1s Y, Is of rank und 5 Ss give for partitioned matrices .1 of order pi 
Satlisty uy 


yrder mal 4 ‘ bye \l lel . 
for some matrix VY of rank and the corre- 
4 ] \ ] VB spondihig partitioned matrix Bs of ordet pl The 
relationship between matrices satisfying (3) and 


) is discussed more fully i! the next section 


where (/ Y)isa, prmatrix which is the dire 
product of the identity trix of order p with the Other reduction formulas for partitioned matrices 
rratrix .Y° ic ire contained in the results of B. Friedman [3], 
5. N Afriat [1], and J. Williamson [10 These 
\ formulas will be shown to be special cases of theorem 

\ » 

The reduction formulas given in section 5 are 
Y \ applied in section 7 to certan special partitioned 
matrices. In particular, a 9&9 matrix, involving 


an arbitrary paramete! re whi h arose in a problem 
on lattices points at the Bureau ol Standards, Is 
reduced by these formulas to a set of matrices of 


order 4. 2. 2. and 1 and from these reduced matrices 


ill erzenvalues are found as functions of 6 








2. Inner- and Outer-Related Partitioned 
Matrices 


The formulas (1), (3), and (5) may be related by 
the matrix Y. For instance, if the matrix A in (1 
is a circulant matrix, 


~ @..-- Gy-2 G, 1) 


@n.-14@-..- @, 


V1) Go. - « Api A ; 


and X is the matrix of its characteristic 
then the related matrix satisfying (3) will have each 
of its submatrices a circulant of order n, and the 
related matrix satisfying (5) will be in circulant 
form; i.e., 








vectors, 








(Ay 1 Loan 
fe—-1 Ao « A, 
A 7 
VA, A, A Z 


A matrix satisfying both (3) and (5) would be one 
which has the form (7), where each submatrix A 
has the form (6). D. E. Rutherford [9] has given 
reduction formulas for matrices of this type as well 
as for matrices related to tridiagonal matrices. 

We will call matrices satisfying (3) “‘inner-related”’ 
partitioned matrices, or i.r.p. matrices, those satis- 
fying (5) “outer-related,” or o.r.p. matrices, and 
matrices satisfying both conditions ‘‘doubly-related,”’ 
or d.r.p. matrices. 

We shall show that dual reduction formulas exist 
for o.r.p. and iL.r.p. matrices and indeed that an o.r.p. 
matrix can be put into i.r.p. form by a simultaneous 
permutation of rows and columns. 

Since the formulas (3) and (5) are stated in terms 
of direct products, the general definition and three 
of the properties of direct products are given in 
section 3. 


3. Direct Products of Matrices 


W. E. Roth [8] derives some interesting results on 
direct product matrices and gives the general defini- 
tion for rectangular matrices, together with some of 
the properties. The three properties needed for our 
results are given after the definition below. Most of 
the properties of direct products can also be found 
in MacDuffee [7], although he deals only with square. 
matrices, 


A, and ra 


ng Watrix 


The direct product of an mXn matrix, 
p&q matrix, B, in that order, is an mp 
defined as follows: 








fa ; a Bb . a B> 

a, B My» ; a B 
ixBb 

Vd, B Ano B ; a BD 


Three properties of direct products 
1. If the products AC and BD exist, 
(A B)(C& D)=ACX BD 


2. If A and B are nonsingular, 
AX B)'=(A"XB 
3. If A is an mn matrix and B is a pq matrix, 
AX B=P(BXA ), where P and () are obtained 
from the identity matrices J, and J] by 


simultaneous permutation of rows and columns 
4. Similarity Between O.R.P. and I.R.P. 
Matrices 


Tueorem 1: Jf A is an outer-related partitioned 


matrix of order pn, urith square blocks of orde) DP, 
satisfying (5 then A can be transformed hy simul 
faneous permutation of rows and columns into a? 
inne r-related partitioned mati ir. A. with blo ks of 
order n. satisfying 3 
Proor. Suppose A satisfies (5 Using property 
3 of direct products, viven above, there exist non 
singular permutation matrices, ? and Q, such that 
AME oF XA S) 
Then, from (8) and (5), 
A[PU,* X)Q|=[PU, * NOQIB, 


P“AP)I,XX T,* X)(QBQ 


So, if we let 


A=fP-‘AP, B=QBQ 


the relationship (3) holds for A and B, and, sinee P 


is a permutation matrix, A is obtained from “1 by 
simultaneous permutation of rows and columns 
From this theorem we can reduce any o.'.p 


the 


and 


matrix, if we know the reduction formula for 
i.r.p. matrix, by first permuting the 
columns and then using the i.r.p. reduction formula 
However, it is usually much simpler to use the 
direct reduction formula for o.r.p. matrices which is 
viven in the next 


rows 


section 





5. A General Reduction Formula 


2 A of order 


satisfying 


whe it A 


have a matria 


orde? pr 


THEOREM Suppose we 


n satisfying (1), or a matrix of 


. > 
for some corres ponding matria B 


} al } 

gsaqvrensy r matrix of rank 7 Then in each case 

there exists a matrix P. de pe nding upon NX, such that 
/) holds, so that B contains r (or pr) of the roots of A 
The only new result here is the formula for a 

matrix sutisiving 5 2S the formula for a matrix 


satisiving | Was vivel by ( oddard and Schneider 
t], and the one for 
by the 
Inatrices are 
the relations nmong them, but the prool is not given 


u matrix satisfying (3) was given 
The three transformation 


howe er, in ordet 


author in [6 


viven below. to show 


in detail as it uppears in the other papers 


Proot simce PS is an wv yr matrix of rank 7 
we will assume 
Pf 
A ( ) q 
where .Y,; is a nonsingular (7 matrix, for if a 


permutation of the rows of VY is nv v to achieve 
this, it can be matehed by a simult i! 
tion of the rows and columns of 41: 1.e., 


NX is equivalent to left 


permuta- 
hutation 


of the rows of plication 


by a nonsingular matrix Q, but a matrix A which 
satisiles | also satisfies 
WAY (A ON) B 10 


matrix 2, 


} columns 


will use the 
as follows: The first 
rmiatrix b which we assume 


In each cuse below we 
whi hy we construct 


ol R ure 


those of the» 


to have the form (9 The remaining n—r columns 
consist of an? / block of zeros and the identity 
matrix of order 7 r, 1.€ 
Ps 0 
o ) ; 
VV / 


Using R we ean the transformation 


matrices for nonpartitioned hiatrices satisfying l 


construct 


or partitioned matrices satisfying (3) or (5) or both 
(ase ] Von partit oned matrices 

Suppose A is an nn matrix satisfying (1) where 
VY is an? matrix of rank ~ and Bis an? j 


It was shown by the author in [6] that } 


Aif we let P R 


matrix 
holds for the matrix 


4 / R r. matrices 
If A is 


blocks ol orde} M. 
the matrix / 


( ase 


2 partitioned matrix of order np, with 


A also satisfies (1 


satisfy Ing |. 


where A replaces the matrix LY, In 


this case it was shown by the author in 6} that, if 


we let Fr / R, then 


P"AP=A 


where A can be put into the reduced form on the 
right side of (4) by a simultaneous permutation of 
rows and columns In this Cust B is ol order pi 
(ase 5 0 R.P matrices 


If A is order np, with 


blocks of ordet P, satisfying 5 
1) with Y replaced by 


a partitioned Inatrix of 


then A also satisfies 


where by property ? of direct produc is, since x is 
nonsingular, .Y, JZ, is nonsingular. Thus, if we let 
Xixl 0) 
P=RxX I ( ) 13) 
A / / I 
1) will hold in this case also, and B will again be 


ol ordet pl 


( ase iD) R.P matrices 
If A Isa matrix of order np with blo« ks ol order 
n. and satisties both (3) and (5), Le 
1(xXXxl XXI,) B 
21h 


where Y is p } and Yisn N A can be reduced 
twice For, by case 3, if we let P R /,. where 


PR, is the p p Matrix given in case } we have 


P;'AP 


where Bis of order nr and JP) is of order 7 o—F). 


It is clear that premultiplication and/or post- 
multiplic ation of a partitioned matrix A, with blocks 
of order n, by a matrix whose blocks are scalar 


multiples of the identity matrix of order n produces 


2 matrix whose blocks are linear combinations of 


the blocks of A 


Also if any set of matrices has the property that 


then 
Ss (ce, A, M=DD (e, X B 
Ol 
Sc, A) X=X(Die, B 
so that anv linear combination of the set satisfies 


an equation of the form (1) with the sane matrix X. 








Hence the matrices B and D satisfy 


BUIXY)=(1XY)E, 
DI,-,-XY)=CU K<VY)G, 


5 
where F is of order rs and G is of order s(p—? 
ry ¢ . . 

rhen, by case 2, there exist matrices P, and P?, such 


that 
Bk * 
P;' BP, ( ) 
0 ig 


G * 
Pr; i » ( 
0 H 


Thus, the roots of A are those of the matrices /, F, 
G, and H, of order rs, s(p—) 
‘(p—r), respectively. 

It would be easier in this case however, to multiply 
the two transformation matrices 2; x J, and J, h,, 
where FR, is the nXn transformation matrix corre- 
sponding to the matrix in case 2, and obtain, by 
property | of direct products, the single transforma- 
tion matrix, P= PR, R,, which could be used to put 
A into the reduced form, 


rh 8), iP and (vi 8 


6 Ff * 
P-'AP 

0 0 G 

00 OO O HY 








For example suppose 


A B 
( dD 


where’ A, B,C, and D are 3> 3_stochastic Inatrices, 


1.¢. 

3 

. 

Sa a, 
$ 

2 b, 
—_— 

j=1 

» ee Cc 
l 

2, a,,=4 i=1,2,3 


then S is an i.r.p. stochastic matrix as discussed in 
reference [6]. 
If we have also 


A+ B=C+D, 
then S is a d.r.p. stochastic matrix. 


10 





Since a stochastic matrix with row-sum s has one 
root s with corresponding eigenvector (1, 1, 








(cf. A. Brauer, [2]), the matrices 2, and 2: would be 
| 0) () 
| 0 
R ( ) R ing @ 
| l 
: @® 7 
Then 
R 0) 
P=R,XR, | ) 
R R 
and 
R 0 
P R R ( 
R R 
So 
R A Is R R BR 
PsP ( 
() R>'(D—B)R 
fa 4 } 7 
() ] B () B 
() 0) d—b 
L 0 () () DD B,D 
where 
a My. a a 
( ) 
ad i. ad ad 
and B, and J, have similar forms 
Another exaniple using Lr.p. and o.r.p. transfor 


mation matrices Is given in section 7 


In each of the above cases, if 7 n, the matrix P is 
much simpler; i.e., for 1, P=, for 2, P=JIp«X 
for 3, P=X™ J, and for 4 if s p P=XxY 


by Williamson, Afriat, and 


Friedman 


6. Theorems 


It was shown by the author in [6] that Williamson’s 
theorem 10] on partitioned matrices follows fron 
theorem 2, part 2, theorem 1 in 6)), 
and the matrices B 


where } } 
ure triangular 


S. N. Afriat [1] gives a general discussion with a 
number of theorems on the determinant and charac 
teristic roots ol partitioned matrices, A | in 


which the submatrices A,, all commute. One of his 


results Is a reneralization of Williamson’s theorem 


and shows that if the roots of the commutative 
matrices A,, are X (k=1, nN; i] l, sp) then 
the roots of the partitioned matrix A A are 
roots of the matrices X NG Since, as Afriat 


proves commutative matrices can be simultaneously 





reduced triangular form by the same unitary 


transformation, it is clear that A A satisfies (3 

where .Y is the unitary matrix in question, and the 
corresponding matrices, &,;, are triangular of order 7 
So Afriat’s result on the characteristic roots of A also 


follows from theorem 2 


to 


Theorem also provides au new prool tor the 


following theorems by Friedman [3] 
THEOREM | Let A,, Ao, , tl, be n-dimensional 
square matrices and BR = b,. BR he p-d mensional 


matrices Suppose that A . A . . A ha ea 
and the corresponding e1ge n- 


A, re spectively Then p oT the 


square 


common eigenvector, ©, 


alues are d,. Xo. 


edger alues of 


( BixXA BX A 14 


l he ¢ ger alues of the matrias 


D=\, B,+rB 2B 


As Friedman points out, this is a generalization ol 


the well-known theorem that the eigenvalues ol 
BR A sare the products Au B ms 9 a 
p where A, are the roots of A and ug; are the 


roots of B 

Tueorem Il. Let C be the 
theorem 1 Suppose the ring ge nerated by the matrices 
Aj, A has an r-dimensional representation } n 
/ hich the re presented by M Ther p? 
of the ¢ igen if the matria 


. : 
matrix considered 41 


matrix A 


genvalues of C' are ¢ alues 


p= 8 Bx M 
theorem | as a special 


prool that 


theorem II contains 


we will only 


since 
with 7 l, theorem 

PRroot Friedman 
statement. “4 
Lo 


prove 
ith his 


Is represented by V/ 


the 
’ 1s equivalent 


shows 


NAM c=] t 


where the columns of the matrix VY are the basis 
vectors tor the subspace 
Let the elements of the matrix P, be b | 
p I t 
Then from (14 (isa partitioned matrix of ordet 


bloc ks ol orde! 


pn wit 


( 2, 0V A gx] 7 
Krom (16 
> 55 1A DOR (XM 
¥ > 62M p IS 
So if we let 
Dij= 2 6%) M, 19 


} 
i 


is the matrix defined by (15) and we 


(1S Lo 


D 
17 


then D 


, 
have, from , and 


(20) 


XD,  (.j=1,..., 
the matrices C and PD satisfy (8 Thus, by 
there exists a matrix P? such that 


dD 
Ok ) 


so that 
) 


part 2 ol theorem 2, 


(21) 


P“CcP=( 


so that pr ol the roots of C are roots of D More- 
over, if the matrix Vis known, )) and £ ean be found 
explicitly by using (11) to construct P? 


7. An Application to a Lattice Point Problem 


The example below Is a matrix whi h nrose inh a 
problem involving lattice points at the National 
Bureau of Standards 


rF4 | 8 7 


3 1 | | 3 
’ 3 l 3 | } l 3 | B 
3 | j 3 l 

| 3 j l 
3 l 3 | ! l 








L B | | 


If we partition into 3 } submatrices we see 


shows that a tridiagonal 


D. E 








matrix, 
& h 5 
h ad h 
ha bh 
/ 
h a b 
v b al 
has characteristic roots which may be written in the 
form A a-+-2b cos 6, where 6 her) 1, A l . 
The characteristic vector corresponding to A is 
sin @,. sin 26 sin né 23 








The transformation for an i.r.p. tridiagonal matrix 
was given by the author in [6] and Rutherford also 
gives all the roots in the d.r.p. case. 

The example (22) falls into none of the above 
categories, however, since B #8"; but we show how 
the transformations discussed in section 4 can be 
used to simplify the matrix so that five of the roots 
may be found without having to solve an equation 
of degree higher than 2. I am grateful to D. E 
Rutherford of the University of St. Andrews for 
giving me the other four roots. 

We first use the transformation matrix, P= (VY x / 
where X is the matrix of vectors (23) for a 3° 
tridiagonal matrix; L.e., 


and J is the identity matrix of order 3. Then we find 


| V2 1 ic" 0 
B=—P;'.AP;, ( 1 C , 
D ( A+7,2T] 
where 
0 Bb UO 
> »7 
( BI] 8 0 B 
2 
0 8 0 


Since A and also the matrices A+ ,2 J are all 
tridiagonal, we now use the transformation matrix 


P,= (1X) and find 
D—-y,y2I ©" 0 
6=P;' BP, ( D CG 
0 C D+4y,21 
where 
P=diag (4—y2, 4, 4+, 2). 


If we write the rows and columns of 7 in the order. 
826435791, we obtain the matrix 


with 


GH a £ 
"or! Mle we 
HG L* M 











where 
4 V2 D 
G ( ) I] / 
' re i—y2 
{ } 3 ) *} 
K B 4 8 |, 3 | op 
L 3 ay Lod 
and M=diag(4+ 2,2, 4—2,2 Thus £ is an o.r.p 
circulant, so the roots of # are those of G+/H/ and 
G—H, which are 4+ 8+y,8+2 
Since A is tridiagonal, let 
4 
rm) 
/ 
Then 
(4 y 20 B/y 2 \2 ) 
t 
P i, t \ 28 ) yz ) V2 
B/y2 B/y2 t 2y2 
. B/y2 B/y2 i—2 2) 
Thus ft is a root ol this matrix, and when we 
delete the second row and column, we are left with 


a 44 matrix. I was unable to reduce this further 
by any relatively simple similarity transformation 
However, D. E. Rutherford the 
four roots of F 


has found othe 
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Bound for the P-Condition Number of Matrices With 
Positive Roots 
Philip J. Davis, Emilie V. Haynsworth, and Marvin Marcus 


s« ptembe I 


Let a matrix A have positive roots 0 
I t P lit unber of he P A 
mmetrie functi of the roots. 
1. Introduction 
Suppose ia a Ss ou nonsingular n ym matrix 
with roots. A» ] } ordered so that 
() \ \ rN 
kor wide class of matrices 1, the ratio A \ 
fives a rough measure of the probable acecuracy ol 


the computation of the inverse of A, or the solution 


of the svstem of equations lor which A is the matrix 
of coefficients 

This measure was evaluated in some detail by 
von Neumann and Goldstine [5|' and has been called 
bv J. Todd [2, 3, 4] the ?- condition number of the 
matrix, Lt 


P 


In oe 


to tlre rect prog al ol / 


neral the accuracy of the results is in proportion 


I> 


I) this paper we shall show that if the roots of A 
nie all positive with 
0< X, <Ad* \ | 


and the determinant of A, det A, is known, togethet 
with one other svmmetric function of the roots, 
upper nd lower bounds for 7? can be given in 
terms of these two know! quantities 
More precisely, if the characteristic polynomial 
for tl atrix c1 Is 
f j 
| ( ( {) n 2 
and we know « and « det A . there eXIst positive 
uppe ind lower bounds for the ratio 
\ 
P : 3 
\ 
terms ot ¢ and ¢ 
To sin plify the notation in the statement and 


prool ot the theorem We use the following device 


» ‘ 


7. 1960 


i << 2 < Upper ind ver bounds 


ven in terms of det A and one other 


hi 
(;.)° 
ol p by This 
equival nt to dividing each element of the matrix by 
‘ We will call the matrix “normalized 
with respect to ke’ on simply the normalized matrix 
the condition number, 7, is the ratio of two 
roots, ? will not be affeeted by such a transformation 
Let D, be the determinant of the normalized matrix 
Then 


and divide each root is 


hew 


Sines 


(i) 


DD 


From Hardy Littlewood Pélva | the nun bers S 


as defined in (4) satisfy the inequalities 
b 


where all the X 


henes 


equality holds only if are equal 


by 5), D,- # 


2. Statement and Proof of Theorem 


Tureorem: Jf we have a matrir A with 


J d a 
characteristic p lyn omial 7) a f the constant term 
! 


LE 


and the Ith ¢ eth ent are he) the follou ng hounds 
hold P. 
,-l+y1-—D : 
Pp. —_ x 
D i ccnall coal 
here D is defined by a 
If i I (ue , the trace i the ma give) we cai 
mpl e the ipper h mnd 
I—D ne 
P- \ D 8 
l \ | D ( 
Proor: Part 1—-l'pper bound. We prove (8 
first, as the method used for this can be applied in 


the prool of (7 


Since P s is equivalent to 


, eq 


tX,A 
| Ds: 
X 


13 








so from (5), we must show that 


(Ay T r,)? 


Dar. S( a) 


n 


(10) 


But the right side of (10) is (A,/@,)" where A 
is the arithmetic mean, and G, the geometric mean, 
of x, hy 0 

Thus, we can prove (10) by induction if we prove 
the following: 

LEMMA: 


(11) 


/\ 


(4: ' (fen) 
a) A\AGant ” 
where A,,, and G,,, are the arithmetic and geo- 
metric means of Ay, . . ., An, Anai- (The A, are not 
assumed to be ordered with respect to size for this 
lemma.) 

Proor oF LemMMA: Simplifying (11) we see we 
must prove 

Ang 148 < Atti. (12) 

But (12) follows from the arithmetic-geometric mean 
inequality (cf., Hardy, Littlewood, Pélya [1]}), 


a,+a -G.41\ ™ 
Gide... Ani <( ); 
a n+-| 
if we let a,=A,, (i=1, . . nm) ANG Ges:= dos. 
Thus (11) holds and (10) follows immediately 


using finite induction on n. 
To prove (7) we note that, by the inequalities (6 
g” I 
Di < Di=— 


X 
, 


so it suffices to show that 


1a 1 — 
tei Lia(!\ on (is 
tr, A, Ly Sy 
Now let w=A,',..., u,—A;' and we note that 
O0<mi< ... <u, and (13) becomes 
(ui + Mn)* 8) (a a) 
a < —— . (14) 
fin, (3: (my ; i) 


But (14) is the same as (10) with the y’s substituted 


for the \’s, hence (13) holds and the general upper 


bound (7) is true for all values of k. 


Proor: Part 2—Lower Bound. From the in- 
equalities (6) 
S, x 
D,=—".> 
Ss; ae 


where, by the definition (4) 


81 <y.8_p > APA, 
Hence 


which proves the lower bound in (7). 


3. Some Remarks 


In reference [6] T. Kato gives an upper bound for 
the P-condition number which in our notation is 
pes. 
D 
This is larger than the bound given in (8 
We derive several from (7), (8 
Suppose that we have a matrix with positive roots, 


({). ives st 


consequences 


which is normalized so that e¢ 


Then ),.—det A, and we can write 
l : , Pp: 1+ \ | (det A ; - 
(det A l—7J] det A)’ 


This implies that 


det A->0 


if and only if P 


det A! if and only if P—] 
Hence, the det A behaves essentially as the reciprocal 
of the P-condition number and may be said to con- 
stitute a reasonable condition number of its own 
for such This lends 
popular feeling that for a properly defined class of 
normalized matrices the smallness of its determinan| 
is accompanied by difficulty in inversion 
A second observation relates to least 
proximations. Express the general problem in the 
language of inner product spaces. We are required 
to solve | y—Zf_,a,r minimum, where y is given 


matrices substance to the 


square ap 


and J, are independent elements. The normal equa 
tions have matrix Fe where (p,q) designates 
the inner product. We can assume zs, normalized 
(hyd r,*=1. The Gram matrix Li yd Is 
positive definite svinmetric and henee falls withi 
the scope of our inequality with k=1, ¢,=n. Hence 
if Gis the Gram determinant ry Al we have 

l <p 1+-41 G 16 

™ l—J1—G 


The quantity @, which acts as a “measure” of linea: 


independence ol J). ia therefore also serves 


as a condition number for the normal equations 
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Some Computational Problems Involving Integral 
Matrices 


Olga Taussky 


Nove 
In this mainly expository article dealing 
tional problem of Gi Pall related to finite 
detail Numerical results obtained on SEA 


1. Introduction 


All computation ona machine with a fixed number 
of digits arithmeti modulo integer and 
therefore ean be regarded aus number theory From 
this fact, in “analytic” problems there arise all the 
difficulties associated with “‘round-off,’’ 1.e neglect 
of insignificant digits. However, there are often 
considerable difficulties in handling actual number 
theoretical problems, e.g., because the numbers may 
Crow large rapidly Codes which permit numbers 
of arbitrary length have been written, but they have 
to be monitored carefully 

Matrices with integral elements have been studied 
number of 
The 
whole vast classical theory ol quadratic forms comes 
under the theory of symmetric matrices; many new 
problems have arisen in vears Some of the 
problems can be divided, perhaps, into four topics: 

1) Inversion a 
2 Eigenvalues and eienvectors 
role in the determination of the class number of 
algebraic number field / by 
it can be shown that the number of id al classes coin 
cides with the AS, 
where A is a matrix with integral elements which is 
a root of the algebraic equation determining /, and 
S runs through all unimodular matrices with integral 
elements QM. Taussky [2, 5 

3) Enumeration of matrices with special propet 


Ls some 


for a long time and an enormous 


problems arise, both theoretical and practical 


very 


recent 


Rosset ] 

These play a 
an 

For 


see ¢ 


rational methods 





number of classes of matrices S 


see 


ties. Here we have the various problems concerning 
latin Squares, the problems concerning matrices 
with 1 as elements and orthogonal rows (see e.g., 
R.E.A.C Paley 1} and J. Hadamard [5 The 


study of firite projective reometries has led to new 
probli His in this field and SOTHEe of them will be dis 
cussed here late \ generalization of these problems 


by the study of block designs which are of 


Is Yivel 
lmportance mn statistical work ser MM Hall OI, 
H. Rvser [7 R. C. Bose [8 W. S. Connor and 
WH Clatworthy 

1) Maximization probl ms. These come actually 


ilso under the problems discussed in but have 


miber 


p 


15 


28. 1960 


irticularly with recent problems a con 


ojective geometries discussed it 


nd SWAC are 


discussed 


been studied rather extensively, particularly in recent 
A prototype among these is the assignment 
problem: given an n matrix determine the 
permutation ?(7) for which 


vears 


ti a 


S ‘a 


— Pi) 


~ 


B. 


Maximum -. Motzkin [11], C 


Tompkins [12], H W. Kuhn and A. W Tucker [13 


sec { 


2. Some Problems 


However, many other problems arise which cannot 
easily be brought under these Of these 
mention a few Matrices with rational integral 
a generalization of rational integers 
and any number theoretical process can he applied 
them any two n matrices of this type 
have a common left divisor, 
common right divisor, a least common left multiple, 
a least common right multiple If Dis the greatest 
common left divisor of A,B, then 


headings 
we 
coefficients are 
Lo 


e.g n 


vreatest a greatest 


D=AP-+ BC 
can be solved where ?P, ‘are agam 7 nm matrices 
with rational integral coefficients Diophantine 


equations can be studied, even the Fermat equation 


can be considered for matrices 
Completely new problems arise as well because of 
the \ quite elemen- 


h 


ol matrices 


noncommulativily 
a 

tary example Is to express a matrix ( , witha, 
( d 


he 1 in terms of the 


b.« ; d rational integers and ad 


l l 0) l 
generators ( ) ana ( 
0 ] ] () 


fic ult problem Is to express when possible a positive 


) \ much more dif- 


definite 7 n matrix C’ of rational integers in the 
form (C'=.\.Y’ where LY is again an nn matrix of 
rational integers and .Y’ is its transpose If the 


& such a decomposition 
Actually find 
Y can be a difficult computational problem (see [7] 


determinant | ¢ l and n 


is known to be possible alwavs to 








3. Hermite Normal Forms 


For many of the problems concerning matrices of 
rational integers the so-called Hermite normal form of 
a matrix A plays an important role. It isa triangular 
matrix obtained by multiplying A with a suitable 
unimodular matrix 1’ of integers and is of the form. 


Chi) hi» h, a 
0 he A>. 
UA 
LO 0 0 hud 








where hy, >0 and 0<hya<h, (ick). 

Another important normal form associated with 
A is the Smith normal form which a diagonal 
matrix. It is obtained by multiplying A with two 
unimodular matrices (’, Vso that 


Is 


fh, 
he 








and h, hy ho. 


4. Application of the Hermite Normal Form 
to Finite Projective Planes 


An important use of the Hermite normal form was 
made in the study of finite projective plane geometries 
(FPP). An FPP isa set of k elements called points 
which are distributed over & lines such that exactly 
n+1 points lie on each line and exactly n+] lines 
pass through a given point. Two distinct points 
are contained in one and only one line; two distinct 
lines contain one and only one point in common. We 
have k=n?+n+1. The number is called the 
order of the geometry. With such a geometry ther 
is associated the so-called kk incidence matrix 


yi 


Is 


A (a 
where 


if the 7th point lies on the jth line, 
if the ith point does not lie on the jth 


16 





Hence any two rows (columns have aone lh Common 


in exactly one column (row) and there are n-+-1 
elements a | in each row (column 
The incidence matrix A is a matrix of integers 


If the lines points ure numbered hh a different orde! 
then A is replaced by PA (AQ) where P and ( are 
permutation matrices. The following two properties 
of A are particularly important 


det A yO Farle (yn l 
AA'=A’'A 


where J is the unit matrix and J the matrix consisting 
of 1’s only 

G. Pall 
forms of A 
free values of n 


the Smith and Hermite normal 
14] He showed that for 


the diagonal elements in the Smit! 


studied 


sec 


square 


form are 
cs | n,n | 
and the diagonal elements of the Hermite normal 
form are 1, or divisors of n, or ( 1) times a divisor 
ol n 
Two different geometries cannot have the san 

Hermite matrix. Thus the Hermite matrix of a 
geometry is rather important G. Pall furthe 
investigated the p-rank of A, Le the maximu 


number of linearly independent rows of A, mod 


If p 
then the p-rank of A 


where pis a prime factor of 7 Is not a tactor 


of n (e if n is square free 


is exactly 


For all other Pp 


p-rank: 


holds As special examples the 2-rank for the 
unique geometry corresponding to n $f was Col 

puted. It is 10. The 3-rank corresponding to 
n=¥Y presented a more difficult task There are 
two different geometries of order 9 The 919] 
incidence matrices of these reomnetries have as 
3-rank 37 and 41 while n?+n l 1) These 
values of the 3-rank were computed on SEAC and 


on SWAC 


The computation on SEAC was carrice 


out by M. Newman; it took 20 hi The Hermit: 
normal forms of the geometries of order 7 5 and 
7] 7 were also computed ; it took 3 hr for this task 


5. Further Facts About the Incidence 


Matrix 
Incidence matrices have been used to study the 
still unsolved problem: For what value of n is there 
a finite projective plane? So far it is known that 
for n=p’, p a prime number, there is always 


finite projective plane. This plane is uniquely 6] M. Hall, Projective planes and related topics (California 


Institute of Techr ology 1954 


determined for 7 ae a ae ey ee This was estab 718. 2. G t \ 
ow *} é ' Lvser, eometr 1 | ( | far m rices, m 
lished iol ‘ by AY Hall 16] and ior 7 ‘s by Mat! Monthly 62. 95-31 1955 . 
computations carried out on SWAC by M Hall 8) R. C. Bose. Mathematical theorv of the svmmetrical 
J LD Swit and R J Walker 17] : factorial design, Sankhva 8, 107-166 (1947 
Kor? 6 there is no KPP whatsoevel see R H o] W.8 ve ~" = ne! A se "A “Mat om — _ 
’ +1 partially balances lesign " Stat 25, 
Bruck and H. J. Ryser [18], 5. Chowla and H. J 100-112 (1954 
Ryser [19 On the other hand for n=9 it is known 10) J. H. Curtiss, editor, Numerical analysis, Proc. 6tl 
that there exist two different geometries; one of these Svmp. Appl. Mat! Ne Yor N.Y., 1956 
11 I S. Motzkin, The assignment problem, Numerical 


Sens Se SeSeen Desargues, the other does not analvsis, Proce. 6th Svmp. App Math.. pp. 109-125 


For } 10 it is not known whether a KPP eXIsts New York, N.Y.. 1956 

{ valuable contribution towards the solution of 21 C. B. Tompkins. Machine attacks on problems whose 
the problem ol determining the possible orders ol a iriabies are permutations Numerical analvsis, Prog 
FPP is given hy the following theorem of Rvser 2) 6th Svmp. App Math., pp. 195-211 (New York, N.Y 


1956 


Let 1 be odd and let k=1 n+I1. Assume nI+J 13] H. W. Kuhn and A 


, W Pucker \ i f Matt Studies 
YA there X is a matrix of integers Then there 24, 28. 38 (Princetor 


N.J.. 1950. 1953. 1956 


erisits a FPP pith j ] points o7 each line 14) Gy Pall, Some theorems , fi te projective planes 
, " . inpublished manuscript, 1953 
| oa on eo = earlier. r odin: rather 15 \ J Hoffman, M Newma I ( Straus and O 
difficult computational problem of demiding whether nn ge oy Tia Mar sects Ataris = age Nap 
a given positive definite matrix of hict ordet split atio Pacifie J. Math. 6, 83-96 (1956 
up In the form VN’ with YY an integral matrix 16] M. Hall, Uniqueness of the projective plane with 57 
| nts, Pros Am. Math. Soc. 4, 912-916 (1953 


17] M. Hall, J. D. Swift, and R. J. Walker, Uniqueness of the 
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Math. Soc. 66, 327-345 (1960 


trices with integer coefficients, J. Research NBS 49, 


1 «) () eorem of stir r and differ 
i von , “29 i ind Macl , 22) H. J. Ryser, Matrices of zeros and ones, Bul. Am. Matl 
Canad. J. Math. 1, 300-302 (1949 P 66. 442-464 (1960 
0. 1 igs ns \ apuce — cm p aeh : 23] R. C. Bose, S. S. Shrikhande, and E. T. Parker, Further 
; , prs — ~ sults on the onstruction of mutually orthogonal 
\ Math, } IS7-19 Ne \ N.Y 1956 - nd the { , . | : AB tae : 
1] R. I 4. C. Pa 0 I g .. M ; gases aa sera * ee poe Fe 
P 38. 311 1) 92 -~ 4 _eeaaaaaaaadl ork), Cana 1. Mat 12. 189-203 
) J H mara Re itlo 1 rie est relat ‘ LUN ; . 
hoon ay B Se. Mat 2) 17. 136-142 (1893 Papel 65B1 13 








Computational Problems Concerning the Hilbert Matrix’ 
John Todd °* 


No nber 21, 1960 
I} " iction | we I il mathematics and pi ictical computational expire 
ent llustrated by i ssl of recent work, at various centers, concerning LLilbe 
i 1X 
l. Introduction If gy To,J so that <r; is convergent 
this gives 
Recent work. In various centers, concerning the pap f mi / wx 
Hilbert matrix shows very well the interaction 
between theoretical mathematics and practical inh ¢ ontradiction with (2.1 
computational experiments Some aspects ol this lhe quesuion as to whether there is an sz, neces- 
work will be discussed here in an expository mnannet sarily not in 2, for which the equation 
Although we shall be concerned mainly with the 
finite segments /7, of the (infinite) Hilbert matrix Hr= ma 
I] i} ) | 0.1.2 is true was raised by O Taussky 20) This 
question will be discussed in sections 5 and 6 below. 
many of the results can be generalized, e.g to 
matrices of the form 3. Applications 
I] a h O,1.2 Among the areas In applied mathematics in which 
the Hilbert matrix or related matrices have turned 
ol up are aerodynamics A. R. Collar [3)}), diffraction 
ol ele tromagnetic waves W Magnus and I Ober- 
H1\0 i 0,1,2 U< 6 | hettinger |12]), and statistics | R. Savage and 
kK. Lukaes [17]). 
2. Theor . 
y 4. Computational Problems— Inversion 


The literature of the last half-century contains 
many theoretical results concerning the Hilbert 
matrix Of these onlv the bast inequality is 
immediately relevant [S p. 254 

2 ow convergent the 

S* S224 ) ) | e >a 
Wi note the following consequence ol this 

y ie The equation Tha rr cannot be satished by 
any ectol ; I] lhert space, / 

Suppose 

Hae=ra 
is satisfied for a vector se/ The it Is permissible 


to multiply both sides by ~ . tuke inne! produc ts, 


to vet 


Tl we'd 
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Theoretically, the problem ot inverting the finite 
matrices H] 
explir it representation of the inverse, a result which 
is essentially, due to ¢ ‘auchy The matrix is of the 
as a double-alternar t, and the inversion 


Is solved because of the existence of nn 


form known 


has been discussed e.g., by A. ¢ Aitken [1, p. 134], 
A. R. Collar [3], and G Polya and G. Szeg6é [15., 
pp. 98-99, 299 300] The explicit form of the 
inverse, and the actual numerical values for 
} 1(1)10, have been given by |] R Savage and 


EK. Lukaes [17] 

Practically, however, t 
esting For instance, although the Hilbert matrix 
itself turns up in idealized situations, it may be 
expected that slightly different matrices turh up in 
the prac tical circumstances which are approximated 
It therefore, desirable to study the 
inversion of the Hilbert matrices //,, 
enjoved a reputation for ill- 
instability,”’ before 


he problem remains innter- 


these Is, 


by 
Tried hanical 
aus 


especially they 


condition 1.e., ‘numerical long 
the advent of high-speed computers and the basic 
papers Neumann and Goldstine [25] and 


Turing 24/, in order to compare the obset ved inverse 


of von 
with the theoretical one 

The Hilbert matrix and related ones 
10] matrix 


Lotkin 


inversion 


e.g... 


have been often used to test 








programs prepared for high-speed computers. All 
the experiments of which we are aware [2, 16, 23] 
and in these several different schemes for inversion 
were used—completely confirm the bad reputation. 
The results quoted by Todd [23] are fully repre- 
sentative. For machines with 10 to 12 decimals, with 
floating point arithmetic to about 8 significant 
figures, there is rapid deterioration of the quality of 
the alleged inverses which are produced——and usually 
complete failure to produce any inverse for //, with 
n about 8. 

A detailed examination of the condition of //7, has 
been given by Todd [23]. Given a matrix A and a 
specific process of inversion, the error in computing 
A™ can be found, in theory. In practice what is 
wanted is a cheap estimate of this, in terms of quan- 
tities associated with the matrix and e.g., the word- 
length of the machine. Various condition-numbers 
have been introduced to measure the condition of 
the matrix. One of these, the P-condition number 
is defined as P(A)='A\/\u) where \ are the 
greatest and least among the absolute values of the 
characteristic values, of A. It has been shown [23], 
that @ being a certain positive constant, 


and lu 


P (71, ) 


O(e@") 
while for the nn matrix 


‘a 


— hw 
t = 
o— 





i—2 1 
4 1—2 


associated with a second-order differential equation, 
we have 
P(C,)=O(n2). 

Using these results in the general error-estimate ob- 
tained by von Neumann and Goldstine [25], for an 
elimination and comparing these theoret- 
ical error-estimates with actual observed estimates, 
shows that the two are in 
(see Todd [22] 


process, 


reasonable agreement 


5. Computation Problems—-Characteristic 
Roots and Vectors 


As soon as programs for the computation of char- 
acteristic roots and vectors were constructed, it was 
natural to test them on the //,, since input pro- 
grams for this had already been constructed. There 
was no reason to suppose that it would be particu- 
larly difficult to handle—since it is considered that 
the presence of multiple roots, or roots close in 
absolute value would be the most likely source of 
trouble. The theory of positive matrices due to 
Perron and Frobenius (see e.g., Wielandt [26]) guar- 


20 





antees the existence of a single dominant character- 
istic root with a characteristic vector all components 
of which are positive, so that the power method 
should be satisfactory. This is indeed the case 
the dominant root, and the corresponding vector 
have been obtained for values of n up to 200 at 
various centers |5, 18 

If more characteristic roots and vectors are needed 
the Gantmakher and Krein [6 theory of completely 
positive matrices guarantees that the characteristl 
roots are all simple and positive and the compo 
nents of the vectors exhibit a characteristic pattern 
the first has all its components positive and each 


succeeding one has one more change of sigt 


and 


The last criteria shows up the weakness of thi 
results obtained For instance, a reasonably effi 
cient using floating point with 8 significant 
figures which gives the dominant root A, of // 
rect to SS gives the smallest root negative! Actually 


code 


Col 


det H.==3.75 K 10 


and 
2.08 10 rv 1.12 10) 
so the results are not too surprising 
The usual method for handling this case is at 
approximate diagonalization of the matrix using the 
Jacobi process oO! the Variation on this clu to 
W. Givens [7] in which the matrix is reduced to 


triple diagonal form and the roots located by i 
So far few actual results of the 


use of these methods have been published 


Sturmian process 


6. More Theory 


Following experimental computations, O. Taussky 


19 showed that, if A denotes the dominant 


characteristic root of J7,, then 
T () | log / 

so that A T The actual approacl to wis ve 
slow, e.g., It is known [5,18] that: A 2.08 \ 
2.11, A; 2.14 nN 2.18 - o } a * 2.27 
the exact order of magnitude of x d, does 1 
uppear to be known 

A constructive attack on the problem alresn 


mentioned in section 2, of the existence of a 
satisfying 


Hv=m 


appeared difficult and an affirmative solution to 


was obtained by T. Kato [9 who used an indirect 
approach. His method, which will now be ce 
scribed, was suggested by obset vations of the mono 
tonic behavior of the coordinates of the dominant 


when these are norn 
Is unity 


vector of f/ 


characteristic vectors of // 
ized so that the first coordinat 
the 
responding to A, Is 


Suppose characteristic 


Numerical evidel cee sugvests that both the sequence 
A and all the sequences : ’ 1.2 


The first fact is known (e.¢ 


ure monotone 
Collatz [4 and the 
more ceneral result 


second Is i special cuse ol il 
regarding the dominant 
a. b of two positive mnatrices A, B 


and # is B re lf 


vectors 
where 1 is aa 
the 


dominates s1 in sense that 


a<6 and if « b,./a is a nondecreasing function 
ol us long us it Is signific ant then 4 dominates 
a in a similar sense, provided further that 2 has all 
its second-order minors positive This is proved 
using the fact (already noted) that in these cireum 
Stances the dominat t vectors can he obtained hy thre 
“power” method The special case used is that 1 
which A4=/H/ / I] 

Now assume that A, is bounded that this is the 
case when we are dealing wit I] follows fron 
2.1 Then lim d \, where A= in the Hilbert 
Cust It follows that the sequences 
ure bounded for comparing the first coordinates 
In the equation 

/] N 
we have 

i Doh \ hn SA 
which gives. for all, 
N(/, 

ry the Hilbert case this Is 
If we write lin it ca the be shown that 
the infinite veetor satisfies 

/] z 

TI e result obtained by Kato 4 while ar swering 
the original question, itself suggests manv more 
For Instance, is the vector the only linearly In 
dependent vectol ( rresponding to \ 7? Are 


the Se ¢ values whi fh OX eed rv? 


haracterists 


7. Recent Developments 


| During the last lew vears the theory deseribed 
i! section 6 has been cle veloped co siderably by 
Tosio Kato and \Iarvin Rosenblum and Sore ol tne 


1 


questions raised there have 
the relevant the following 

M. Rosenblum [28] has shown that every complex 
number with 
root of //(@ 


by using tl 


hee inswered Among 


results are 
a positive real part is a Characteristh 
for 6 fixed, @>0 
ie theory of spe inl functions to obtain a 
vector explicitly In (29) Rosenblun 
ed the spectrum ol H1(0 fol 
anv real 6#0 l, y 4 this was accomplished 
the Titehmarsh-Kodatra 
differential operators 


This was ace omplished 


characteristi 


completely determi 


theory of singulat 


Meme 
tisiti 


I Kato [30] shows that if 17(@) is the Hilbert 
bound ol H/ o), 1.e \V/ “ Tr COSCE rt 0) i « I 
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V/ 


teristic value of HH (6 


vec 
wh 


H 


(6 is a 
with a positive characteristic 
tol value A \7(@) 


ich has a positive characteristic vector 


0 T @>4, then everv A charac- 


and there is no characterist 


2. A reference to Hilbert’s own work about 
Is 27 where he evaluates det I] Kor sore 
related results see R. B. Smit! 32], W. W. Sawver 


Oo 
trot 


The problem of matrix inversion 1s discussed 


n the experimental point of view in [31], which 
udes further results on the Hilbert matrix 
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Index to the Distributions of Mathematical Statistics 
Frank A. Haight’ 


Feb 
\ f ‘ rhe j refers 
bef ] ir 1958 is pr ‘ ! nate 
, a) ¢ j 
t ‘ I i ‘ 
li ! covered ar it La 
listribu l 
he is Divar i in ct ( I 
ot } i! ! 
lox 1 we 


I. Introduction 
1. Background 


The author began this index in April 1954 with the 
Auck 


limited intention of supplving his students at 


land University College, New Zealand, with a 
small reference pamphlet It appeared that no 
textbook of mathematical statistics contamed a 
complete treatment of all the distributions that a 


Thus aun index Wis needed 
In 


and 


student might encounte! 
to facilitate selection of the upproprhite book 
1955 the collected results were mimeographed 
during the next three years several hundred copies 
of this early to 
requests 

Then, at the invitation of the National 
Standards, the author spent the summer of 1958 at 
the Bureau’s Statistical Engineering Laboratory 
supplementing and editing thre index lol publication 

lorm 

course of this work, 


version were sent out in response 


Sureau of 


In the present 

In the 
SOUTCEeS } uve been lhe luded and the journal coverage 
extended through 1957 Information 
supplied by readers of the original version has been 
mistakes have been 


a number of additional 


bye he a 


incorporated ana various 


cCorres* ted 
2. Organization 
Thy, 


material given under each distribution consists 


of a number of entries. most of which are prov ided 
with one or more references In the case ol the nol 
mal distribution with mean m and variance \ No 


number of entries ts fairly large, and there 


standard ordet 


1.1 


lore 


the 
the is most easily seen 
| Functions and parameters 
Ly. Derived distributions 

a) of linear quantities 


h) of quadratic quantities 


15, 1960 
ol statistical distribution | ) ad 
rial given for each distribution is a list of er 
s which characterize the distributior lerived 
potheses, (e) miscellane he 
mit II], binomial, discrete, distribution 
tions over miscellaneous 
! IV} imber of entries \ I 
lred for the normal dist 
j fy Its alread th 


(". Estimation 
a point 
b) interval 
}). Testing statistical hypotheses 
a) by linear statistics 
b) by quadratic statistics 


I: Miscellaneous 


These cutegories are by no means always used for 
important distributions. With the limited in- 
formation available, a complete listing of the head- 


leas 


ings in such cases would be wasteful as the majority 
would be empty. Keeping in mind the above order, 
it should not be difficult to find the required entry. 

entry will be indented; such an 


as a continuation of the pre- 


Occasionally titi 


enury should be read 


( eding one 
3. References 


The references to the literature are of the following 


types 
A. Coded 
a) Journals, e.g., |c}4:17, which refers to page 
17 of the 4th volume of the journal designated as 
c] in table 1 


h Books. e.f.. |? 53. wl it h relers to page 53 
the book designated as [12] in table 2 
b | neoded, O2.. Trans Am Math 


contorming to the usual volume and page relerence 


ol 
Soc. 17:382, 
stvle 
(". Reviews 
a Vathematical Reviews is designated by 
h) Zentralblatt fiir Mathematik is 
by Z 
MIR and Z references will in no case offer a of 
a paper appearing in coded journals and therefore 
mav be considered to indicate publications in obscure 
from the poml of view of the present work) sources 
\loreovet every effort has been made to avoid MR 
or Z references to an uncoded paper quoted and very 
few duplications of this sort should be found 


MR, 


designated 


review 
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The choice between direct (i.e., coded or uncoded 
and indirect (i.e., review) references is frequently 
available. The one that was actually inspected is 
given. All direct references were collected before 
the search of MR and Z. Consequently, each entry 
corresponding to a direct reference is based on the 
paper, and never its review, and each entry cor- 
responding to an indirect reference is based on the 
review and never the paper. 

As it is difficult to distinguish priority in a large 
number of references, a chronological table of the 
coded and review references is provided. This table 
also shows what volumes have been systematically 
searched in the preparation of this index. 


4. Criteria for Inclusion 


A. Distributions. As a general principle, a dis- 
tribution is included if its density (or probability 
function is a known, explicit function. The follow- 
ing exceptions may be noted: 

(a) Certain families of distributions are men- 
tioned, e.g., Pearson and Koopman, whose densities 
are specified only implicitly. 

(6) Certain distributions are mentioned in 
terms of their cumulative probability function or 
characteristic function. 

B. Entries. The general principle governing the 
selection of entries is that each entry must exhibit a 
property of the distribution in question. Excep- 
tions to this rule generally take one of the following 
forms: 

(a) Historical information about well-known 
distributions, although not systematically sought, 
may in some circumstances be included. 

(6) Important applications, such as those which 
led to the discovery of the distribution, are usually 


supplied. 
(c) Bibliographies. 
Reference to tables has been excluded in almost 
every case. 
It is clear that applications must be severely 


limited. With a slight exaggeration, several whole 
branches of statistics may be considered applications 
of some particular distribution, as exhibited for 
example in the following table: 

Distribution Application 
Quality control 
Analysis of variance 
Probit analysis 
Random processes 
Applied mathematics 


Binomial 
Normal 
Lognormal 
Poisson 
Deterministic 


5. Relationship Between Distributions 


In some cases (such as 2.1 and 2.3) the relationship 
between two distributions is asserted in their desig- 
nation. In others (such as 5.3 and 5.15) a very close 
connection not pointed out. In most cases, 
however, koown relationships are simply listed 
among the miscellaneous properties of both distri- 


Is 


24 





butions. In choosing between these procedures, an 
attempt has been made to follow current statistical 
usage and terminology 

Very similar principles have been used to decide 
for or against independent listing. If one distribu 
tion is relatively important and its equivalent much 
less so (for example Chi-square and Erlang), they 
are listed together. In 
reputation seems to justify independent categories 

It must certainly be supposed that many of the 
trivial distributions of section 8 could be included in 
some larger category or even combined with each 


other cases independent 


other. Such a systematic classification which would 
exhibit all even if worth 
certainly removed from the purpose of this index, 
and has hardly been attempted. 

For example, it is well known that no. 8.1 con- 
tains as special cases all the distributions of sections 
| and 2; very likely it also contains dozens of others 
listed. Nevertheless, to indicate this situation by a 
system of subheadings, applied to all entries, would 
quickly undermine the utility of the whole work, 
since it is the special ease, rather than the wel ral 


connections, doing, is 


principle, that occurs in statistical practice 
6. Notation and Terminology 


In univariate distributions the stochastic variable 
is always denoted by x, in bivariate by x and y and 
in multivariate by x,, Xx, quite regardless 
of the domain of definition. This 
from that of certain authors in two respects: 

(a) The letter n is not used for a discrete variabl 

(6) The statistic obeying a particular distribution 
For example in Student’s 


departs 


usage 


is not used in the density. 
“t” distribution, we write 


rather than 


This practice is justified not only by the need for 
uniformity but by the belief that the alternative i 
wasteful of the alphabet: t, F, z, D, Similarly 
we prefer to call distributions by the names of them 


rather than by 


discoverers (or reputed discoverers 
the symbol used to denote some statistic found to 
satisfy them. Of course, all known designations are 
found in the final index 

In many books the expression [(x 
denote a probability density However | 
monly used in mathematics for an arbitrary lune- 


Is employe ad to 


is Col 


tion D(x) has therefore been selected as a more 
distinctive means ol representing this sper ial 
function, 


In the discrete case D(x) replaces the probability 
distribution, which is often written p,. C(x the 
cumulative function. 

When we come to the characteristic function the 
situation is a little more complicated Using t ior 


works venerally have to 


characteristh 


the variable, statistical 
define several svmbols fo 
various quantities, tor example 


functions of 


x (ft characteristic function of distribution of x 

» (t characteristic function of distribution of 
nx 

E(t characteristic function of distribution of 
Xx. etc 


As we will be dealing with manv different statistics 
and possibly their characteristic functions, it is more 
economical and systematic simply to abbreviate by 


the following system: Ch(x), Ch (nx), Ch (X), ete 
Thus it is not necessary to select a new letter to 
denote the characteristic function of each new 
Statistic. 
However. this practice leads to equations like 
(Ch X c 

which sities be offensive to some, howevel cle ul the 
meaning Such readers are advised to interpret the 
equality sign as an abbreviation for the verb “‘is.”’ 


This inte rpretation has anothet portant connec 
with the being used \ variety of 
been emploved to deseribe the 
variable distribution. 


tion notation 
verbs has 
between a 
for example 

x obeys the normal distribution 


relation 
stochastic and its 
wit! mean m 
and variance v 
x follows the normal distribution with mean m 
and Variance \ 
x 7s a normal variable with mean m and vari- 
anee Vy 
It seems equally felicitous to assert this relation 
ship by the convenient abbreviation 


D(x N (m,N 


which may, uf preferred, be regarded not as a mathe 


matical equation but as shorthand In anv ease it 


makes possible an unambiguous condensation of 
the facts 

Similar remarks apply to the expressions MGF(x), 
FD p which are used to mean moment renerating 


function of the distribution of x and fiducin! dis 
tribution of the parameter p 

(nother application of this use of the ec unlit sigh 
relates to the svmbols C R p MILE p MME Pp), 
UMVUE(p), BANE(p), and is exemplified by the 


following 
(’.—R.(¢ v/2n 
MLE( s 
UMVUE(r) =ns?/(n-1 


For the meanimnes of these and other abbreviations. 
the reader is referred to the following list 
D(x Density or probability function of a 


stochastic variable x 


C(x 
Chix 
\IGF X 
PGF X 
KD p 
1h) 
x 
\ 0 
Mx 
Ky 
I 
p 
GM 
HM 
MN 
; 
b 
3,,6 
( R p 
My7E(p 
BANE(p 
LR 


seg 


OC 


Cumulative distribution function of x 


Characteristic function of distribution 


ol x 


Moment 
tribution of x 


vyenerating tunetion Of dis- 


Probability generating function of x 


Fiducial distribution of a parameter p 
Mean of a population 
Mean of a sample 


Variance of a population 


Variance of a Samy le 
( Ss ee X X 1 


kth central moment of a population 


kth moment about the origin for a 


population 
kth cumulant 
Correlation coefficient in a sample 
Correlation coefficient in a population 
Median 
Geometric mean 
Harmonic meat 
Number of items in a sample 
Slope of regression line in a population 
Slope ol regression line In a sample 
Asymptotic large sample 
Pearson's betas 


lor vari- 


param- 


Cramér-Rao lower bound 


f estimates of the 


nnces ¢ 


evel p 


Mlaximum likelihood estimate of the 


parameter p 


Minimax estimate of the parameter p 


Minimum Chi square estimate of the 
parameter p 

Uniformly Variance un- 

biased estimate of the parameter p 


PHvdTLaUdn 


Best asvmptotically normal estimate 


of the parameter p 
Likelihood ratio 
The likelihood function TI D(x 
Sequential 


Operating characteristl 








BCR—Best critical region 
(Q—A quadratic form 
K—Expectation 
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necessary 


II. Distributions 
1. Normal 
1.1. Normal (m, v 


A. Functions and Parameters 


D(x)=(1/V¥2av)e"""" v. 6\108, 5\34 Al? 
(S]97, (9243. 

Ch(x)=exp(— }vt?+ mit 1]211, |5)\62. 

MGF (x) =exp(4vt?+tm): (6]172. 

Derivatives ete. [\d|2:181. 


Transformations: \c]39:290. 


Obtained from Pearson’s differential equation 
[/1¥o 
14) 4 Fe 

Called Type VIL: 11\45. 

Limit of binomial: [4]58. 

Variance of X and s: [3]42. 

Var(m;)=6v*n™', var(my)=96v'n7! and many other 
constants: [2|224. 

Calculation of constants and numerical examples 
[11]88. 

Mean deviation E\x—m (2v/xr A9T8S80 
(1]258. 

Probable error = .6745¢ |4]58 

am, = (2k—1)v*: S5|NXTIT, |8\98 

Quasi-range : [\d|24:603. 

B. Derived Distributions 

D(x)=N(m, vn='): [9|270, |6|10.2, |2|243, |4]100, 
\w] 1:93. 

D(x/s): [3] 139. 


D{(n—1)4s~'(x—m)]=Student (n—1 


(21239, [4]112, |w}1:74. 


6)217,|5]98, 
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D{s~'(n—1)!X], where x; are samples from N(m 
m, not all ‘qual \d|19:406 
D|(x—m)/range]: [d|22:469. 
D (range) ete.: MR13:762. 
D(>k.x,)=N (Sk ym, Dk?v, 4199, (892 
x—m/o)|? is Chi-square, (x/o)* is nonecentral Chi- 


is Type II], product of normal 


lol 


In} 
Is Bessel, quotient ~normal 
18|1—150. 

X,/Xz 


| N(0,1),  m.>o [2]253 


0,1 


square, 
variables 


Xx 
lar cr¢ 


lisa 


m,—m, 


when m m., and y \ \ 





[6]263 


D : 

n)s 
Student 
[c]29:350, [ 


pf (" . ay ny (Xx 
Snedecos - 


Dik™ 


) smallest sample value 


“1 


n hh ‘ 
, contidene e ellipse 


1|374, 
No 


value from top d\25 


S40 g42:408 


‘ 1 
Gamma (4 ] 


VUR4:164 


D3 x—a 
D(HM 


( ‘oothi ent ol 


a\95 OH 
13] At), 
Typ V. 


Dirange/s A\17 :366, \ 


, ‘precision constant” mome 


d\3 At). a\9? 132 


etc. 


D (logs 


1) ns-\ Xx 
For unequal v, 


D{n,(n.—1)vos7/n 
4\115, 1O\197. 


n l)y Snedeco! 


— 


Testing and confidence intervals 


Power function 





For equal \ ( R.(o v(2n)-'. hence s is efficient 204 
D(n,siv Nos3v x7(n,+n.—2), Efficiency of estimates of « 137-182 
n,(n | In estimating o 


[70|1 is yon 1 . \ “= X m Ss more ethic ent than 


Diy = (n L)s X AW 16. 


D(variance ratio k\11:136. vin—Tany C84 
Distribution of various statistics from k normal 
populations with common variance: e\17 22 ( R a/m): e)8: 04 
Ranking variances gl51:621. n/(n—1)]s* unbiased for v W152 
Distribution of many quantities in a wide variety of Estimation of o for industrial quality control 
cases J. Soc. Stat. Paris 96:262 ee ie 
= a x MEI 14a Estin nation ot ¢ trom percentile S Lt) gr 
D(x,s 2238, MRS:161 Estimation of v and a A\18:584. |bl1-78 
Db.) for n=4 is hvpergeometry c|250:411 Cc R. (vy 9v’n-|. hence s? efficient 484 
Dtd,) for n=4 is hypergeometei ¢]25:207, | MLE(v)=s 6110.3 
D(midrat 121-100 MLE(s depends on whether 1 know! 3134 
LD rane Cll7 oF el{X-] 2 / 12) Estimation oly el12:5%7 10:85 
Quasi-range D28:179 In estimating v. ~efficienev of s 1 
+4 b\15 193. If m constant, s* not sufficient 177. [45138 
1) NS N26:11 In est utling Vv, ~efficiency of mean deviation S76 


1) XX af 7 A\IS:265. Unbiased estimates {1 M.S. Tranelat n No. 98 


D : he J #70 (‘losest estimates of m.\ 33:45 
KI) I \ X. VI y 407] 14. \p 231 | unbiased estimation ol ment absolute adeviation 
> 7 ~e” 
x 3:56; 
FDIinis~! (x Student (1 on VR13:3 
Tests rLYUS 0133:173 90) A9:279 In estimating m,\ ~varlance-covariance matrix 


snaty OP mr HA d\10:68 X, Ss are oments estimates ot ! \ ] iS 


FDC Helmet 1, 1/sy'n 189 My Rr Ye ait 


Ranki PANS ; 3] X s oimt sufhierent 12% ell? :211 
IS STS] 


R if | 
- Estu ation ot m d\T3 . 110-321 el8:12 
Xx / Ustimation 135 :1S¢ 
( KK VI x eflicient j V2 13 Estimation ( ] 
, 40) 17 LS 2 Nf Wu x estumation , 17 ) 


were as hited MLE from censored sample 22124 


UMVUE(n x SOT. |t)4:167 ( ensored sample 0139:260. \e\43:225. MR14:569. 
- } 1° VRS P41 


x unbiased pli :1a0 
: Estimation of m 1145:214, [d\17: 43 
Minin ix estimates of 1 ad ?] 718. j d? IS I 
; If v known, confidence regions for m are x+k,v! n= 
} 9? se 4et . = 

Best CeNSItY unbiased estimate of m A290 29S where k, are the p% values ol tiie normal lol fe 
vefhiciency of & is 2/r=.6366 1|490, |uj22:706. Confidence intervals for m (150, |6)224, |LOTSI 
Estimation of m when it must. be integral. ete seq confidence intervals for m [d|18:427, [b|19: 


Meat of k \ ilues irom top and bottom has ~effi- (‘ontidence intervals for \ 
cweney Zero 1490) Interval estimation of vY and o el6:117. 
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Confidence limits for v,/v2: [4]131. 
Confidence limits for m,;—m, with same v: [4]130. 
Interval estimation of m,;—m, (Behren’s Problem): 


[3]91, {d]18:601, [d\14:35, [d|15:430, [d\20:616, {d\21: 


507, \e\4:39, 108, [d\24:390 {p]7:232. 


Confidence limits on m and s: [e]2:13, [o]8:83. 
Confidence intervals for m, v: (6]227, [3|79 
Tolerance limits: [d]13:398, [d|17:208 

(~) [d]17:238, [d]27:171, [t)1:164. 
D. Testing 
Testing m, Mo: [dj14:149. 
Tests on m: \6)259. 
Unbiased regions for testing m,=m,: (3320. 
Power function of m>mpo [3305 
Test of m using range in place of s: [d\17 371. 


Hypotheses on m: [1}533, [4] 149. 

Tests based on (rectangular a priori) distributions 
of m and v: Z18:158. 

Seq. tests on m: [d)16:171, [e]10:364, 
[c]37 334. 

Control charts on m: 
Testing m, against m, 
counting methods: 
Seq. hypotheses on m: 

Student’s hypothesis: 


368, [b|9:250, 


(6) 16:131. 
and go» against 
[e]17:80. 
|\d]20:502, [g|40:308. 
[g|31:318. 


a, by quick 


Student’s is best for testing m,=mz: [3]285, 291, 
[e]12:79. 
LR test of m=mp is Student: [4]150. 


Power function for Student test: \d]17°192. 

Seq. Student test: [c]37 :326. 

Comparison of two means: 
[e|35 ‘88, [o]4 331. 

Comparison of k means: [t)7:1. 

Testing whether variance is constant: y\20:114. 

Three decision seq. test of m: [y]23:22, Amster- 


dam Mathematical Centre Stat. De p. Rep SP34. 


[10]190, [e\88:252 


Testing whether many means are all zero \e]8: 
70, \4)176. 

Testing X;— X»: [2]29:21, [c]41:361. 

Testing X,—X, without assuming v,=v, [d|9:201, 


[13]433. 
Linear hypotheses: 
Joint tests: [t]6:25, 
Testing outlying observations: 
H: o=a: [3]287, |\d]8:193. 
OC for x’ test of o= ao: [d]17:179. 
Seq. test on a: [e]10:369. 

Seq. test on o,= a2: [e]12:63, [b]11:101. 
Tests on v: [6|267. 

Fisher is best for testing v,= vo: 
Testing homogeneity of variances: 
The most powerful test of 


c]27 :161, [3|292, 300. 
So. 


| 
7 


[e] 17:67. 


[3|289. 


(c]31 Py 250. 


H: oo, m 


Alt: o=o,, m=m, 

Hypothesis of equality of many normal variances 
(16:89, (c]29:124. 

Significance of smallest of set of variances: 


117. 


s|10: 
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Or 
3278. 


Critical regions for m and \v 


Bibliography of testing equality of variances 
a|109:457. 
Seq. ratio test termimates [e]8:342 


OC function g\47:191. 
Power functions of tests 


[d]17:189. 


Whether two samples are from the Same norina 
population [n]7:3, [k]17:302 

Tests for normality [e]28:295, [e]27:310, 333, 
(c]34:209, [o]1:125, [i)20:152, Z19:74 

Tests of various composite hypotheses [19 A905, 


l¢ 19:30, l¢ \10: 79. 


Decision problems: 


KE Misee llaneous 


Independence of x and s (Student-Fisher Theorem 
[7]19:108, J. Math. Soe. Jap. 1:111, MR1:34¢ 

s*, xX independent: [4]108, MRI4: 

Normality if and only if xX and 


(b]15:55 


he ies 


‘ie 


e2 


independent 


d\13:91, \d\16:400, {d\13:91, NBS Rep. 2267, J 
Math. Soe. Jap. 1:111 
Normality if and only if D(x, s)=L. sn d\14 


197. 


Normality if D(x) Div) =¢yx7+y MR10:125. 
Various characterizations of normality fe] 132359, 
128 :126, 39:59, (|d|27:858, fel 4 180. An J 
Math. 61:726, Math. Z. 41:405, Z13:214, 3rd 
Berkeley Symp. 2:195, MR16:1034 
Generalization of Student-Fisher theorem 2U 
248 
Independence of quadratic forms: [o]1:83, Pro 
Roy. Soe. Edin. 60:40, |u\80:178, [e\37 293, [e|14:1.95 
1€|15:427, |d|20:119 VWR12:509 
Baves’ theorem: ln|16—1:113 
Cochran's theorem: (4|107,68, 
More generally: (\d/|11:100 
History of normal: le|16:402 
History of distribution of s c|23:416 
Distributions which converge to normality: d 
10:247. . 
Discrete analog [el44:365. 
Regression of x and t, where m=a—be™' A\18 
596 
Comparing percentage points of two normals 
ld 19:98. 
k samples: [3|295. 
Truncated sample: (d|23:237 
Sampling from N(Z aym,, v 1\160 
Max 2x;, min =x: ¢\40:35 
See also: [125 :389, [d\1:151, [d\6:197, |d\7:77, 
[d|\10:365, |d\13:235, |d\17:483, |d|20:123, |d\21:362 
1d\21:557, |d\22:596, |d\23:43, |d\23:384, |d\23:547, 
1d|25:16, (N3:309, (e\9:6, |13\3845, |e|10:522, |c|13:287, 
[cl16:239, \|e\24:184, |al79:4545, [1|7:23, |ce\81:238 ( 
34:61, 98, |nlb:3, [n|12-3:65, [e|40:116, [c|32:226, 301, 
\g|48:550, [d|25:6386, G98, |v\5:337, [g\b1:88, Brit 
Assoc. Math Tables (8rd Kd vl Pp rrviui, NBS 
R p. 2545, C.R. Acad. Sci. Paris 238:444, [0)1:88, 
App. Sci. Research (Netherlands) Sect. A 3:297, J. 
Franklin Tnst. PHO ZZOF, N RB S Rep PIC? . Am m 
Math. 57:821, Ann. Math. 35:312, Nat. Acad. Si 
28:297, [y|24:2, 218:225, MR16:52, MR14:1098, 


Z5:A73, Z19:317, MR17 753, |w\7 2193 





1.2. Normal: N (0, v 


D(x 27 exp(— 3 Xx°s L0\50, {e\3 121 
Chix exp I vt AV 
MGF(x)=expd vt N53 
C(x C125 :375 
an = o7*(2k)!/2*k! L0\5 
2nd cumulant =v, others zero [2|67 
Pearson type 21141 
D(x N (0, v/n 21175, |n|10-3:90 
D x/V Cauchy IS\] 150, w|1 27 
D(x Type II] 6 :212 
Dis AN2AG, [ul28:456, Z15:118 
D 2 Xi/z2V Gamma 4 a 
D(> x ete c\35 2:47 
Dis x DQ, Q mils, 
D Q) \ yl! ) ol rank i AU cigenvalues all 
L:Id9:4e8. lel24: 122 19:357 |e|\380:407 
Ch(Q)=[1l (1 —2itvk 
D{x n? (Xx Type II (29:13 
x) assuming v is Tvpe III d\28:510 
Vins Student testing €\37 265 
Dirange) for n=3, unbiased critical region 3 
FD(\ T'vpe Viz x?/(n—2), 3} (n—2 pl? 22 
C R. (v 2na'y 1\484, |[p\7:159 
If Din No. 3.5 VWR14:391 
MLE(\ ee \141 
s* UMVUE of v, but s not of o 152,54. 
Neyman-Pearson on hypothesis testing c\20: 
178 
- H ‘5 
\lost powerful test ol Is x s*.< 
L Alt J 
> ff 
Completeness edo! 
Unbiased critical regions 3|212 
Testing s7/s5 etc e\5:157 


Testing 31:106 

Various devices for showing urea 

Inference A\15252 

As **Maxwell-Bolzmann” 

Variance of Gini’s mean d 
2\217 

Mean difference C\Z8 2432 

Properties of f(x), where x is N(0,\ e|17:211 


serial correlation 


distribution 12\39 
llerence 1s 


see aulso: / 12 255 Clo oll, J Ai 178. ( >] 
A) ni13—1:51, |u\80:330 VRS364 WR3 2, |aidso 
[24 
1.3. Normal: N(m,]1 

D(x 2n)~) exp[—3(x—m 
K(x?7+x Mor ar(\7+-X 

2— 3a Kx 2/% K(e* 

( Varte ‘ e* Si720 
D(x 3|2 
D(x"), D(x?+-x3), D(x7+-x})! 8}96 
D X\ Bessel VWRI10:200. 
Fk Dim N (X, n 13|85 


distribution (m 


¥ sufficient 3)8, [e|17:211, | p\7:161 


Baves 


X is consistent 3135 26 
< is MLE 1|140, |e|33:125 


X IS minimax 


Efficiency of ‘ 637 >If 
Confidence intervals for m 363,70, | p]7:222 
MILE for x° test A\25 :580 
Remarks on testing 1113:62 
\lost powerful test of H n“~m 
Alt: m=m,m 
Is X >< 3127 
UMP test of m=0: 13\452 
Seq. test cl43:452 
Completeness edo:sls 
Unbiased critical regions 311 
Testing equality of several means e|8:69 
Testing; called ‘Laplace”’ PoD5 I 
Peculiar composite hypothesis on m O06 


Inference h\15:52 


) > , 
Pitman’s method 1924 
> , . 
Rang C1ID 2405 
See also c\27 66 C191 2202 €\36:400, ly 7:96 


3rd Berkele y Sym posvum 1:14.97 


Gaussian 


1.4. Normal: N (0, 1 


D(x 2n)! exp 4x 7\129, 
General « x pose Acta Math 77:1. 

C(x) as continued fraction 2)150 
(*(x) as a series eW1OI:16 

sounds on C(x C142 7265 


Property of C(x Vath. Zeit 41:405 


C(x VWR10:267 

~C(x VWUR16:628 

MGF (x) =exp Gt 

Ch(x)=exp tt S|1¢ 1\100 

OO 2k)!/2*k!: d\5:32, (d\11:358, [h\1:13,198, 
11208 


Absolute moments Z1 226 


Median and quartiles Cle } 
Dink Normal: c|19:227 

Dis elb:138 

D(x d\1:340, |e\8:138 

D( 2x x" (n 21231, '4\103, |9\331 
D(>Sk,x d\13:17 

C (ns '§ Student: SSE 
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D(Q), D(Q,/Q,): 
D[=(x,— &)*] 


C(r) expressed as an integral: 


(e]17:37. 

(9333. 

[F339 
[d|19:375., 


x*(n—1): 


D(x;x,) by Mellin transformation: 
D(x,/X») [d\19:375. 
D(,s): 


D(range) for n=3,~D(range): 
[c]5:318, |o]8:155, [c]36:142. 


Moments of sample median: 


[c}35:120. 


Cauchy: 
(ol? :65. 
[e|3 aed, 
[d|26:600. 
D(extreme deviate): 


C(range) : [c]32:341. 


D(2x3/Zv?) =Snedecor: [d|19:378. 
Ch fens of estimates of v: [d|19:257. 
Generating functions: Z2:200. 


le 36:96. 


[e|3: 


Estimation of dispersion: 

Estimation of mean deviation: 3254, [e|35:3 

Variance of median: (c]23:361. 

Testing N(0,1) against various alternatives: 
[e]30:139. 

Multivariate analysis: 


[3|INXVITI 


Limit of binomial: (71134. 


Central limit theorem: (7127 139, [1|29:206 
K™ value from the top: [1|374. 
Censored samples: [e]41:230. 

Ordered samples: le|11:23. 

Stratified sampling: [d|5:138. 

Variance in two samples: [nJ13—3:49 
Ratio of two ranges: fd\21:112 
Tetrachloric functions: [e]14:157. 
Approximations: [d|17 363. 
Sheppard’s tables: [e)2:17 4. 

Grouping: [7|32:135 

Moments of order statistics: [d|41:200. 


Occasionally called Laplace-Gauss, or even Laplace: 
Acta Math. 77:1, C. R. Acad. Sei. Paris 232:1999 


See also: [d|4:109, [d\17:350, [d\22:425, {d\24: 
133, [d\24:297, [13|63, [e]18:3.95, |c|24:98, jc|24:280, 
[c]25:195, [i)6:209, [17|No. 41, [m\6:120, |d\22:418, 
[yj24:22, [ul]29:231, 24:66, Z28:266, 220:39, 144, 
MR12:191, Z18:412, MR17:756, C. R. Acad. Sei 


sss 


Paris 238 34 tty Phil. Trans. Roy. Soc. London A237: 
231, MR7:18, Z5:366, [y\4:189. 


1.5. Truncated Normal 


C(x): [1|248, [6)243. 


>» 
Introduction, estimation, examples: [15] 144. 
D(2x)): [b|8:223. 


Fitting: [c]3.9:252. 
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Estimating m and v: [g]44:418, [g]472457, [f\9:489 
[0|3:37, MR7:461 [e]40:52 
Distribution of estimate of o 15\316 
Censored sample: [c]42:516 
MLE: (7|32:119 
See also: [d|9:66, |[d|20:458, [g|47:379, Brit 
Assoc. Math. Tables (38rd Ed.) vo 1 p. errv, MR2:231 
1.6. Generalized Normal (Kapetyn 
2rv)~! exp [—4v~' (f(x) —m)*]df(x 1|.93 
C.—R.(m v/n, C.—R. (v 2v?/n 139 
MLE(o In-'S(f(x, m 
MLE(m) =n7'!2f(x, 
See also: le 5D:168 
1.7. Normals Added 
D(x l+-k 2m) t exp h(x-+-m cv '(2a 
exp —}3(X—m))V 1\2:340, |h’\2:63 
Var (x +k l-+m?+-k(v+m 
Method for partition with example 8|5247 
Semi-invariants: (alt7 21 
D(x (7|11:219 
Three normals added d\i -237 
More generally : M31, |d\5:230, le\3:85, MTR 
LS5. 
Sampling theory: n\8—3:67 
Many normals added [e|37 2429 
Called “compound normal”: JI2180 
Bivariates NS:328 
See also ¢)40:460, [e\14:369, MIR11:258, MIR1? 
L102 
1.8. Lognormal (a, m, v 
D(x x—a 2rv) texp[—4y log (x—a)—m 
parameters and moments: [/4|/60, [4|121, [do 
|] I58. le Ay BLY) 
Graphical determination of parameters w)\9:102 
Mean=a-+e' ’ Tare’ "(te l 
Another form 
|  e x—a 
D(x exp log 
2 e(x P 2¢ b 
\- a . 
m—be! 2, mode=be~' a, GM=¢ 
a b, moments, tables, regression, 
examples, bibliography (14:30, |b\7 2145 
Moments, transformations fw)? :152, [w|8:83 
Complete treatment with bibliography: Aitch 
son. J. and Brown. J. A. C.. ** The Loqnor mal [is- 
tribution” Cambridge, 1957 VWRIS:957 
Estimation of m e\10:341, 
MLE: \g|46:206, Intl. Congr. Math (1950) 1:481, 


Regression ANF 3196 

~ Tests on m A\28:1044, 

Called “Galton-Macalister’’ CWIKk Lc 0d 

Called Gibrat: Kendall and Buckland 
ary of Stat. terms.”’ 

Used to 
12:448 

Deduced from hypothesis about errors, ete 
28:1 41 


‘A diction- 


approximate Fisher distribution d 


x—a b—xX lorcnormal 17 |No.46 


Transformation c\B6:155., 
Versus Geochimi 


Acta Q-53 


normal (Cosmochimica 


Discrete lognormal c\37 :362 

Compared with normal by means ol Galton-Kape 
tvn apparatus s14:129 

Truneated lognormal 82150, |e\38241 

Lognormal (0, 0, 1), E(x eh v= ‘ 81120 
17 ¢ 17 \No ». le) 22:109. |d\4:30 
See also N14:120. '1)138:161. lel12:121 hie 
/ hit1:19 1175:182. \e 179. lel22:1 46 ai3 
762. aif As (2:57. le\86:1454 cl38:427 a \ 
G00. J. Frankl I / ) a | 50-339. PAOLA 
P51: ID. 2P51:617 mo0-90 ( 3 °/f) ait {9:15 


and Ena Chen {} 


j 

ié 
> } 
,9. Ind 


) W113:29 


Soe A)216:309, J. Phys. Chen 0: d 
m Thy ene 325 J10: 173 VRS Nature 
1.9. Wrapped-up Normal 
1) X kde Bull Sor Vath France fit 


France / 4 
55 :335.1d\18:589. Hand 


Spi hive , 7 7 


32. 67 :1.C. R. Soe. Math 
Sup rT lt 


Berlis 


1.10. Gram-Charlier 

Two-term D(x 2a | —k/6(3x—x*)| exp(— 4x 
317083. 137 

General Gram-Cl irlier Ava) 

Dix Ae HZ 299 

D(x 5 3212 

Dis 653123 

t-cest OH 21 

D(x f(x) N(m,v): d\23:467, [g\26(P) :233 

D(various statistics q\4cl 

Log Gram-Charlier: 28:14 

Type B  Gram-Charliet AS82183 18:57 
Trans. Amer. Math. Soc. 67 2206. |\d\20:376. »*)? 

MGF factorial moments 75:218., 
See also al88:576. |a\89:129. (el383:126. lel3e 
?} C\38:b8 Q9 ( B94 D5. N71 47. / IZ -IS 


T_A.M.S. 67:206. Z18:320, Z22:2/ 


1.11 Bivariate Norma! n(™ p) 
Ir v 
D(x,% 220, 02(1—p’*)'] exp’ —4[v,v.(1—p 
X— M),)*V2—2poe,o2.(X—m,)(vV—m vV—m),)*V 
61165. 15\89.14\60 
Introduction, properties, examples 15\585. 
Another form (Koopman-Darmois e|8:322 
Chix, exp) 1(m,s+ Mot L(v,s?+2pa,o.st+ vot 
1|\287 
\IGF \6 167. 
D(xX,.x Normal bivariate JU] 
1) X v cC\Z: 37 j 
1) V/xX eZ L285 
Dixv) and D(x \ 18\1—151, Am. Math Vonthly 
jf 
Chixy ai42Z:S2 
Diu 120 
If p=0, Da B(1, n—2 10|160, 178 
If p=0, Dib B(1, n—1 10\180 
D (correlation ratio No. 5.3 LOTS] 
~Da c\10:507 
D (rs)Se, $?,83): uj29 26 
Distribution of various statistics e\l 2] 
MILE(m,,m»,v),V2,pe,0 X, Y, Si,S3,rs,s 3 | 37 
Var(x v,/n, var(s v,/2n, var(1 I l—p 
COY Ss) 5 p°o\.o 2 1 COoVvil Ss po l p 2n., 
COV X.v po;do Nn ate! 
MILE from fragmentary information A\3:163 


X, V are joint efficient in estimating m, and my); 
—- = lis ns 1 
x, 1 rss have joint efficiency 
1 lon l on l ; 
I lin ] hi 
x—m,)/o V—m,)/o./and|(x—m,)/e 
V¥—m,)/o,] are independent and N(0,2(1+p 
N (0.21 po respectivels 6\217 
Estimation: d\17 23.95, \e|8:322 
Estimation, testing 1 \60F Gh 7dad 
Censored samples rjiors 
Distribution ratio standard deviations rier4s 
Confidence limits for 1 c\29:157, J. Nat. Inst 
Personnel Res 6:155 
Confidence limits for m,/m d\13:440, MRI13: 
My)? 
Sufficient statistics e\17:212 
Comparison of two correlations 10\203 
Tests of sevel hypotheses on the parameters 
d\11:410 
Testing equality of two r’s A\12:279 
some tests uw 7:46. 
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Testing equality of variances: [e]1:13, (bibliog- 
raphy), \a\109:462. 
[c]25:102. 


Sequential tests of p: 


r—T2: 
[w]8:202. 


Truncation: MR2:231. 
Fisher’s original work on r and p: [n]/—4:1. 
Sufficient conditions for normal bivariate: [e]6: 


3899, MR15:805. 

[c]27:145, 227. 

x/o,+y/o. and x/o,—y/o. are independent normal 
variables: [e|3129. 

If p=O0 then (v,/v2)/[(v,/v2) + (s?/s3)] is Beta: 

Called Bravais distribution: [2]19:3 

Z15:310, MR1:246. 


k samples: 


[e}31:10, 


Properties: 


See also: [d|4:196, (d\14:141, [e]7:260, (g\26:129, 
[c]22:1, (c|25:356, [c]25:392, [f\8:328, [n)9-3:90, |e\39: 
238, [i|27:221, |6)218, [v]5:311, [d\27:1074, |e|44:289, 
(74:85, Harvard Ed. Rer. 1946, p. 52, MR4:280, 
MRS8:283, MR14:1102, MR7:212, |y\1(.No. 4) 20. 


1.12. Bivariate Normal N ( °. Vi p) 


» Vo, 
D(x,v) =(24M!)— exp[— 3M-!(x*°v,—2pe,o.xy + y’v;)] 
where M=|"! 972 v,v.(l1—p?); Ch(x,y) 
po\o, V2 ; 
exp(— $n~'(v,s?+ 2pa,o.st + t’v2) (9|308, [2] 
22, 334, [4)60, [8)116, [15|588, [10)95, 106, |c\80:8. 
Qo=3vi, ay=3pejie2, axn—(1+2p")vi v2: 2|80, 
|4]60. 
Moments [o]3:2, |e]12:177 
Central moments: h’\4:73. 
Incomplete moments: (e]13:401, {e|40:22. 
Product-moments: [ce] 12:86. 
As limit of binomial: |a]91:548. 
If p=0, o,=o2 called “circular normal,” D(r 


(c]29:137. 


C(x,v) with other properties: 


properties : 


Cumulants in terms of moments: 2189. 
Var(r)=n7'(1—p*)?: [2|336. 
Var(b)=n7'v,~'v2(1—p?): [2|337 

Marginal and conditional distributions: [4]62. 
Regression : [3] 144. 

Correlation and regression with generalization 


(7]24:1. 
Contour ellipses: [s]/ :208. 
[o]1:103, |d|18:565. 


(d]14:195. 


Bilinear forms: 
Quadratic forms: 
Ch(x?,y*): [2]336. 
D(variance-covariance ) 


Wishart: 


(2]340, [1]29.6. 
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3 /9 


D(r) 21342 
Dir) for n } 


u| IG 536. 


D(b)=Type VII: 1402. 
Dis S» No. 8.3 if @ 02: e|2:65, c|31:9 
Di(e*.e w)2:155. 
When, further, p=0 generalized Student Cc} 
303190 
Simple function of x/y is normal a)93 :442 
Div/x ¢]32:16, |1]20:61. 
Dixv), Dix/s MR3:171 
Joint distribution of Pearson betas c)7 :38t 
Y(si, S35, rs So) ete: d\5:283. 
D(s?, s €)5:139, |e]25:126. 
~f/ (nh ])3 b D)a 
I ( ) Student (n—1 
1 ( b*v 
) 
D( : o—Me ) Student (n—2 5.13, [2] 
\\ bv, - 
. 4S ; ] yf 
Moments of dist. of covariance from 
fo, 1 
n(?! ,): — |gi8s4 
0,1,! _—" 
Di(radial standard deviation D[n-S(x Xx p>» 
(y v)}}: A15:75. 
D sample covariance 3|359, |10\138 
MLE p): 3|33 
Confidence intervals for p aS1. 
MLE(\ »V2,p 2|339 
Estimation of p when \ V2 A9:149 
Estimation of p by rank correlation d\7 :40, 


cl40:419 

d\21:272 

3]138. 
e]5:151. 


Truncation 

and estimation VRIG:498 
Testing v,/Vv» 
Testing 0;/02 and p 


Test whether two samples are from same popula 


tion: 3)140. 
Hotelling’s generalized T applied to tolerance 
limits d\14:90. 


Odd fact for p=0 d\18:442 


See also 0|12:90, c|2:369, C\4:498, |c]17:176, 
€)20:295, g\27 :254, a|83:128, c\29:7 A n|2 G40. 
[n]7:6, [c]32:196, [c]38:371, |1]7:220, |1]24:17, |n]13-1: 
21, G5, |ul28:457, Z212:267. 

1.13. Trivariate Normal 


D(x): 
Moments: 


10\255 ,260 


ol414, c\40:23 


4 


Correlation 174:158 
Partial correlation C1091 
Yielding 2 2 ? table }()-2) 2 
Student test for partial correlation ra 
Snedecor test for multiple correlatior aL 
See also \-] 12:9 / / 


D(x ( a \ exp Sd(a 
x 1) \ where A f 
invers¢ of corre] oO trix / N20 
(2 
cas ( \ 
C h(x exp 1=m,t LSS p,,0,0,1 it 


\loments CV4ZU LLU VRS LZ 


Marginal distributions, conditional distributions, 
regression AVAL 

Independence ol quadratic forms \ 

Distributions ol moments partial wma iltiple 
correlations : ISS I? :23 IS 220) 


Diproduct moment Wishart C120:5.2 ? 
271, MR10:387, |b|17:82 


AIS 35°33 p> Af) 
Cou ila ts 0 orn i ceneral ed Varlance 
Te) / 
, , ; : : 
Independence: 0 adistt utLO Ol Canis ale second 
order moments 
1) multiple correlation die: TU¢ 
Fiducial distribution ‘os | 


Multipl and partial correlation ( 
Lineal regression theory ) 
1) () x AVEEP/ 
VRI3:142 
Divector correlation C|28 2353 

Sampling AVN, |d|6:202 

MLE(n m X1, X 6\18? 
Hotelling s 


1) Varlous () 


LR test for variances equal and correlations zero 
ayy): 20 

LR test [or independence ol variables dill 17 

(‘haracterization e\l4cot 

Hypothesis ol equality ol means A123 

lh cle pe nadenes oO! sets o} variables 4 

Probability that all n variables are positive d\Z¢ 
\ 

Tolerance regions 27 217 

Testing Varinnee-covariance homogeneity Col 
If. jelo4:311 


Student test f 234, 
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Truncation lj:17 
Various tests 
Quadrivariats C\43 220 
sibliograp! vy of tests of hypothesis 
Varlinces called ‘Bipolar distribution eVl4c6l 
Variables separated into two sets ( 


Many Sil ples CWT 2221 


Qo? ‘ 25 8. olf 9 Ais 0. ( L3 I} 2 
WRI? 
Z10 Alin 


ly \V/ ith Soe 


P-4Y hle-TS7 ao2:200 
78, MR12:345, MR15:141 
15:220, MIR10:312, 7 


2. Type III Distributions 


2 Type III (p, q 

) 
1) \ | \ ‘ ()} 

rig 

\ 
Dix : gamma 4 

by! a 
D(x Type LIT (k/a, k Erlang 4 339 


Chix | itp Ne 1172 

MIGE (x if O\l / 

Y {| p x p qq | y ) (| l oe | Z ?| 
| \ i ly l \ i ly | 

! cumulant =q( 1)! p f 


1/2a 
21400 


“Gp °, He = 8G (Ga =)p 
Vath 


=o =qp "ye 
Arithmetic, geometric means Student 13 


("(x ¢125:379 

Di(x/(x+\ Beta 14\41 

~C(x VR13:553 

Transformations vy=x", v=¢ AO:17¢ 

Normalizing transform Pro Pak. Stat. Assoc 
) j At} 

Transformation \ x-+k 14:118 

Moments, Ch(x), cumulants when x=v—<¢ 18) 
/ 136, 1—144 

Type IIL (p, p+1 e\d:17¢ 

Dix Ty pe III np, nq ANZA cll S:335. w| 
] 7 , 

1) nq l)n x Reciproc al Type I] p,nq l 

~DQAx 135:297 














D(x/yv)=Beta of second kind (p=q): [wl :74. 

FD(p-')=Type V (nx/nq—1, nq—1): [3]87, |e] 
30:408. 

Bayes D(p~')= Rectangular: (391. 

D(xy) where x is Type III and y is Type V 
[27]39:64. 

D(x—y)=No. 8.92, where D(x)=Type III (A, N) 

and D(y)=Type IIL (yu, 1 

D(HM): MR4:164. 

D(xy), Ch(xy): MR16:377. 

D(xy)=Bessel, D(xyz): Studies Presented to 


Richard von Mises, p. 301, Acad. Press, 1954. 
MLE (p)= Moment estimate (p) 


q/r, correcting for 





bias = (nq — 1)/nX, sufficient, ~ efficient, not efficient 
|3)26. 

Var((ng—1)/nx)=p*(nq—2)"!: [1]505, |ul45:214 

MLE(1/p)=x/q: 3]21. 

Var(x/q)=p~’n-'q™"', sufficient: [3]21. 

Sufficient statistics for p: [e]17:212,219. 

Ordered LSE is MLE for 1/p: [d]25:3165 

MLE (p,q), variance-covariance of estimates b) 
14:187, |c]42:22, [r]1:18 

UMVUE(1/p)=X/q, with Var=p~“’n-'q 353 

There is a sufficient estimate of a 3)26 

Estimation: [e]8:324. 

Minimax: [16]64, [e]14:57. 

Gauging: (e]15:192 

Closest estimate: [uj33:217. 

Testing mn such populations: [3]325. 

Slippage tests for p: Koninkl. Nederl. Akad. (A) 
59 :329. 

Testing equality of 1/p: [c]31:205. 

Confidence intervals for 1/p: 74, [e]6:118, 


MR5:128. 
Truncated distribution: \g|4:411. 
Estimation from truncated Type III: 
[d|\27: (98. 


ld PE B59, 

Truncated samples: [e]40:52. 

D (gap between events) generalized Poisson 0 
3:123, [o]11:101. 

Characterized by independence of sum and quotient: 
[d|26:319. 

Discrete Type III: lel44x 

Mills ratio: [d]24:309. 

Normal limit : Am. Math. Monthly 50:98 

Renewal theory: ([d|11:448. 


Trivariate: [d\21:550. 


See also: — [d]7295, [d]8:17, [d]24:407. [e]2:150, 
(b]11:101, |cl21:263, [d|25:640, |c]24:300, |v|2:330. 
[q)50:904, [7|39:171, lg\7 295, Am Math. Monthly 


50:98, MR17 3756. 
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2.2. Type III (p, 1 
D(x) =pe™, “negative exponential’ [53 
Type X. 
a=1/p: [d|59, [2|48 
v=p [5 |67. 
Moments 211.42, [8|100, [10\18 
Characterization: 't17:60, 3d Berkeley Symp. 2: 
195. 
Cumulants: [2|87, [10|40 
C(x): [e|25 :379. 
Mean difference: IC]28 2452 
MGF(x)=(1—t/p 10\37 
\MIGF log x [oe]? 2296 
Grouping corrections: CISI2433 
A priori distributions of p A he 
D(x+-yv): [S|95 
If Dix pe’*, Diy Pe 
Ae x >\ 
then Dix—y 
Ae! X \ 
where A(p+ P)=pP 
Examples and applications 8129, 79,83, [j|20:36 
lel39:168 
D xX : icl39:168 
Dit No. 8.9, MGF(¢ pl? :153 
>x, where x, is Type III(p,, 1 VRS+42 
Rank variates: (e124:210,271, [r)4:153 
Estimation: [e]35:187, [g)48:493, [s]10:167 s 
15, |p|7:152. 
Censored sample: ¢141:230, [d|23:237 
Moments of D(s): ¢|22:538. 
Variance of mean deviation ~4/3np 21217 
Testing p= po: AI:8 4 
Sequential test: ¢141:252, |d|27:460 
Testing against four other possible dist: ( 
Confidence intervals: rIS4. 
Estimation from truncated exponential 7 


Relation with Poisson, D. L. Gerlough and A. Sehuhl 
Poisson and Traffic, The ENQO Foundation. 1Q55 
[o]2:13 


See also: \A\7219, [d\25 2555 IS: 488. [q|50:90 
; | 
Ifx V c, Mean C+» Val skewness= 2, kurto 
Pp p 
sis=6, Ch(x), cumulants: [18]1—-136, 1-144 
MLE: k\8:52 





lo 


LR 


identi al, ( 


test hypothesis that n such populations ure 


t< i) 


Zl oe 


estimates ol 


k samples 


Re st lin ( 


| . 
of >, |Xt4 x 


al mand @ P90 


Life testing: 125 :373 
Original Nevman-Pearson paper on hypothesis test 
Ing C\20:221 
Censored samples 52:58 
Estimation by order statistics [26:58 
Quasi-range A\28:17 
LR tests d\12:301 
seq. testing oli 
Confidence intervals VRS 243 
( H: p=po, c= 
sf] 
1 Alt: p=p,, c= J 
See also d\25:409, |d\2 8, [c]80:402,416 
le 2x9 
Varwell-Bol2man) 
x? is Type II (p, 32 1\3.9,60, |s|1:30 
Connection with N(O, \ 12\40 
D(X n|10-3:90 
See also n\17—1:125 
2.3. Type Ill (1, q 
Dix | e~*x?-', “Gamma” 10|149 
Pq 
r Ciq+r)/Tiq LO|150,163 
k qu 1)! N96 153 
Chix i—1t)™ 17|No. 3 
Skewness q 1O\161 
HM q | 1O\168 
D(nx Type Il (1 ng C19:228, (ni lO—s791 
D(x c|24:47 
D(x,—x P1252, |c|24:293 
D(x, +-x Type IIT (1,qi+q 10\151,|p\7:101 
D(GM) as a series 2|251 
With generalization d\5 2277 
Dx, /(x, +x Bi3q,, 34 10\153 
D(x,/x Beta of second kind: 10|158, 160, |p|?: 
102 


576069 61 4 


D(x,;,—x.) involves a Bessel function if x, and x) are 
from two separate distributions and D(x,/x)) is 
Fisher a7 2:51 

~ DIX N25 -63¢ 

D(x, GMX D(x) D(GM/ X ¢|30:287 

C R q n(ad dq log iq 

X is moments estimate of q, not sufficient 

Closest estimate: 592216 

MLE(q lor GM pir :lt 

K(x q, var (x) =q/n, efficiency (Xx 0) 1|504,4 

~ Dilog GM Normal 1\505 

Confidence intervals p\7 2:2 

Mellin transform 19:37 

Log log (1/x 8:71 

\ B log x C6168 

Multivariate generalization 2 

Tetrachoric functions: e|14:161 

Fermi-Dirac, x—e is Type IT] ) 2.const 12 

) S:701. |s|1:30 
See also A)25 2401, (10\161. 1d\22:425,. |e|24:281, 
0127 2409. [e\80:415. le\86:165. |\d\25:784. le\10:314. lq 
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01 :467, [e\44:265, [d|22:418, Can. J. Math. 3:140 
2.4. Type III (1, 1 
D(x ( ‘exponential’: 6\217 
( hy X | | it 17 Vo 33 
D(X,s 1\3:128, |d\4:133, [d\4:1389,142 
x V c, Type X, confidence intervals 6117290 
D( =x Type II] by convolutions d5:13 
D(>x,/i D(max x, b)14:43 
Doubly truncated, Ch(x n|10:3, |17|No. 32 
Ratio of two ranges: ld 21:112 
+ R. theorem false for x= ( 1\480 3|47 
sor also d 22 2425 1136:152 mib: TZU, [d|22: 
Is 
2.5. Chi-Square (k 
D(x 1 /2!*1°( 4k) |x*®¥“'e“ 96, (10)164, [2|17, 
102, 18\134, MR8:161, |1\3:353, |18|1—161 
Type III(4, 3k 
a k, \ 2k 11|28 105 
a. k k 2 k 2s 2 1\25 
My 2k. uu SK ba {Xk 12k uu 32k »k 12): 
7) At)? 
Cumulants [e]31:216 
Chix 1—2it)-* 8 
C(x), relation with Poisson I\3:357, |c\37:313 








Introduction, properties, examples: (7: 

Obtained as dist. of normal variance: 
148, [pl7:98. 

D(x,-4 

D(log x): 


Xo): [e]7:27. 

[e}34:170. 

2) if D(x) 
1). 


MRI B 15 


D(24xy) =Chi square(2n chi-square(n 
and D(v)=chi-square (n 


D(x,s) for k=2,3,4, n=3,4: 


D(x,/x.) = Beta of second kind (4k,, 4k, 10\177 
D{[x,/(x; +x»)]=Beta (k,, k2): [10\177 
Obtained as D(> x?): [10|169. 
~ D[(x—k) (2k) N{k, (2k ~ D(2x 
N[(2k)!, 1]: [1|251. 
~ D(—2 log LR): [d\9:60. 
Reproductive property: [4]105, |10|177. 
Normal approximation: [(d|17:216, [d\27 :786 


Large parameter: (¢]43:92, Proce. A.M.S. 6th Symp 

> | t I 
in Appl. Math.., p. 251 

Convolution for a pair of Chi-square variables S 
134. 

Percentage points: [e]41:313, [1333168 

~ Significance levels: (d]14:57,93. 

Elderton’s tables: le|1 7155. 

Minimax estimation [16] 17. 

Connection betweer X° test and xX? distribution 


[e]19:215. 


Original Neyman-Pearson paper on hypothesis test- 


ing: [¢]20:263 


Queueing: [b]16:82 


Approximation for small sample: A\9:158 


See also |\d| 18:89, (g|29:37 2. c)22:298, ajss 
87,95, |a|87: [c]29:389, |c]31:346, |e|34: 


42, [c]29:133, 
| 


-/ 
} 
368, \c]32:268, (c]40:421, |p]7:98. 


2.6. Non-Central Chi-Square 


_ eres 
D(x) =e7""e-#2-1 5S " - €136220 
j=0 '(gn-+))2*)! 
' exp{kit/(1—2it) 
Ch(x) pore. 
(1—2it)”* 
k,=2""'(r—1)!(n+rk),~form: 18]1-162 
Logarithmic noncentral Chi-square 0]7:57 
See also: (2]34:57, [e]41:4538, [1]36 (supplement 
18, [d|\28:678 [cl44:528. 
2.7. Helmert (p, q 


(x/q)®-*e-# /* 


qi'(}p)2'?-? 


[5}94 
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For 2q x, p=2, called Ravleigh 12|39 
D(x*/q*?) = Chi-square 
Called seminormal 20:61 
Refs, Chix 17|No. 42 
p— | - 
a p 5 , = q,V pq a 
Non-central d\23:467 
See also c|23 L418. kele« Eeng., Nov L95L D 
L004 VRS:161. J Appl Phys 23:1 37 
t 
2.8. Reciprocal Type III (p, q 
a 
1) Xx | xX ( pq i) 
iq] 
i 
Mi il} p, Val ’ hvpe V, 
q | 
a pq) T(q—1 1) (q+) 2180, 142 
Obtained from Pearson equation 11) 
Various constants with nn exan pl ‘ 
A3:20 
(h(x) in spec lal case 17 |.No / 
If gq :, Pp 1, Chix eXp \ 
s|2O8:318 171 No 49 
Nore rene! illy a7 225 
MLE 1:19 
As distribution of precision constant im normals 
ple a\9? ] 32 
( (x C125 :579 
ser 1isoO C12 IAS CWO 2165 ‘ ) 
2.9. Generalized Type III 
D(x CU(1+x/a)* with moments of D(s 
¢]22:52 
Dix [pr gq? (p)(q+x)Pr ‘ with semi-inval 
iants, D(x)ete ([¢]21:28? le]24:293 
Ch(x ( Ix/p 17|No. 8 
D(x 1 /pP'(q)|[(x—e)/p]*"'exp|(e—x)/p],e<x 
MLE ¢.p,q 3|39 


Variance of estimates |42 
Tables a1:191. 
Estimation q\48:336 


Dix [e~™*b?/T p)je"’*(x-+-a from Pe arson’s equa 
tion [4]74, [2]124, [11]65, 

D(x ln] 10—3:91 

D(x A(x—c)*"'e*-", x>c, p>0, q>0 
[1]249 





Various constants, with an example 1 1\6¢ kp q mode « kp p 6157 


One root of quadratic in Pearson equation ts Chix pe i—p)’ N62, |2]45,108 
| ’ t /) ' 
I1\44 MIGE (x q pe . NAS. (10\88 
Dix X2) W here each is Generalized Ty pr Ill, or one PGF q+ pt)* 18]1—146 
Is Generalized ype 1] and the other N i, 
NI52 ~('(x rte 99 
Baves’ Theorem nilo—lill4 K(1)x g\49 169 
( ounting radioactive particles G1S:260 Dinx Binomial nk,p NPs 
Two Generalized Type 1] distributions added DD 5 €|44:262 


] 


nis Pe ai FD (p 13? 117 


D(x—v) in terms of Legendre functions: AJR14:456¢ 


et ilso d\1:1a0 adii:19] di}z:18s c)1:293 
el3-31] 715:173. |e113:13 cl16:1] nit 2 - OR e132 Reproductive property by convolutions P21 
24 40 2284 \ convolution with respect to p wi6 2168 
| | 
( K p pq kn IAI AN] pli 160 
2.10. Wishart Univariate X is sufficient pl? :162 
kt 2x Is efficient and 1! binses 114 
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| 101153 L4\42 


3d Moments 241! 


Mode p ; pq 1). ete n IS 

a B p+r, q B p, q H1117 bh) 4:1.2¢ 14\4 : 

HNI p—Z p+q-—2 L166 

G\M—exp (0/d}p)|log l4p)—log }(p+q 

C(x A) 20:451 

Moments |?/] 

Obtained as a Pearsor Typ 417 

From at example O45 

C(x 11921, |e] 252379. |e|38:42 ?:208 | 38 
19 

D(z 57 \ 5 

Wher p=q cj19 A530) 

Special cases and variants 2 

Beta k, k 1409 

Beta (2n—2, 4) connected with tolerance limits 
Hoel, Intro. to Math. Stat 

Correlation ratio in samples from uncorrelated 
bivariate normal is Beta (k—1, n—k A352 


1414 0) 


DiI 


c]21:1, 


Beta 


c| PW) 1O\IS1, 184 


q p 


X 


If x is Beta (p,q), vy is Type III 
1, p+q) then xy ts Type III (1,p 


Standardized Beta variable is N(0,1) as 
p—o,q-—a 11|252 


Dixy): — [e}9:36 


L10O\156. 


Roy. Soe, 
154:133. 


y 
P’roe 


A 


No sufficient statistic for (4p—1 
] 


d ] 1:37 > 
h SJ of 


Mellin transform 
nore wv nerally 


log log —; aresin yx 8:71 
Generalized €\43:237, |c|44:441 
Transforn C136 :165 


Connection with Fisher and Snedecor A419 [c] 
21-350. 
Nonecentral Beta d\26:648 
Connection with binomial c|41:30 nilS:12] 
Fitting straight lines y|24:2 
\pproximated by Poisson 0) 1 
Range of rectangular (0.a) is 
D(x n(n—1)a~®*x®~?(a—x 193, [c]25:417 
In trivariate normal analysis 41 
GM (x x,) end GM (1—x , 1—x 
are joint sufficient \4 
In rank correlation AIS 
See also adileé: iS. d 17 » 20 h 1:21 r C\Z2 S84 
C122:391, |c\23:143, lel2@7 74145 10\154, |c\80:140, [e 
> 2] ? ‘ C14 HS. CWO TG ( g yl 27 ] ( 67166, 
15:75, [el39:204, [c]37:219, [c\40:281, |q\48:831, 
|] 
5.4. Type II 
Dix x—m)?*a aB(4,p+1 
a<x<m-—a A141 
Properties VWR10:1 >] 
Transform to Student C)ZS20U8 
Dix nilo—3:9] 
te X cg 12, CIZ) 2) / 
~ Dist. of rank correlation 
coefficient (Pitman €\30:259 
Called Thompson’s distribution, relation with Stu- 
dent’s, normal approximation d\27 :784. 


fume tion: 


m location, a sealing, p shape, 
dA\3:86, 


( R Tables \ 4 


-p , 
If m=0 
Distribution of Spearman’s , 


11 \62 


like lihood 


1h a 


for large n: 2|401. 
A numerical example 
Estimation of center {4°55 


Krom Pearson system 








5.5 Partial Correlation 





r[3(n—k+1) ny 
D(x) , i i(n—k)] | (4 —x2)4m k-2), 
—1l<x<l: (1412, [(]24:198. 
Type II with m=0, a=1, p=4}(n—k—2 
Ch(x): [17|No. 19. 
If k=2, transform to Student: [5\99, 


If corresponding population is zero: 


[10|256. 


parameter 


As No. 5.1 with r=0: Cowles Commission Papers, 
New Series No. 10. 
If population is non-normal: [36:16 
See also: [d|18:81, [n|2:684, lo 3 45 la IP -5SO. 
5.6. Parabolic 
D(x) =[3(a?—x’)|/4a°, —a<x<a, v=a’/5, grouping 
corrections: [c]39:432. 
Estimation: [d|26:505, [m)\6:120 
5.7. Type IX 
m-+ 1 
D(x)= (1+x/a)™, —a<x<0: 
u 
2\142, [17\No. 16 
See also: [d|7 :26, |c|24:234,240,263. 


Type VIII for negative m: (17|No. 23, MR4:21 
5.8. Type XII 
(1+x/p)™(1 x/q)-™ 
D(x) = 
\P/4 (p+q)BU+m, 1—m) 
m l, PIS g¢: F711 43 


See also: [d|7:27, [17|No. 24 


5.9. Correlation Determinant 


D r[}(n—1)]*!xh-® 
As rik YT [4(n—2) ‘Tl3(n—k 
a, = (n—1)'-*(n—2)(n—3). . . . . (n—k), 
Vv k(k lyn -OCn |/ A117, | 4 120 


Downton calls this “Geometric,” and mentions, in 


‘onnection with LSE, the following special cases: 


ce 
[d|25 :304. 
I. D(x) =px®"!, 0<x<l 
C(x) =x, a, =p/(p—1), v=p(p+2)"'(p+1 
I. D(x)=pb-°(x+a)’"', —a<x<b—a 
Il. D(x) =pv~'b~*[(x—m)/v!+-al? 
m—av!<x<m-+(b—a)v!, a,=m, 
a=p!(p+2)!, b=p-'(p+1 p> 1 
log log (x [|8:70. 
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5.10. Triangular 


D(x l 1—x', D(X) from rectangular when n=2 
D(range [c]25:417, [y]24:22, MR 3:171 
Dx 
D(x) »>ksxsk l 147, S132 
2k+1 
Stratified sampling [c]13:48 
x—m 
Dix Vo | + 22} right triangular 
o 
a 2,20 xX m \ 2a d A Pit id 25 :308 
D(x) =4R"(4R— x—m x—m!S4R, best linea 
estimate of m and oa 1d|25:318 
Testing [125 26.95 
See also d\2:48, |d|\28:179 
l 2k 
Dix 1+k »OSxSa 
a a(a—x)_ 
1<k<1. ealled “linear’’ d\4:2 
l+x, l<ox<0 iain 
Dix » called “Tine / 3 
l—x, Os x< 1 
Chix 2/t | cos t 17| No. 14 
) 
D(x a—x), called “semi-triangular 
a 
mn n/.. V a 1S, grouping corrections CISI 2452 
Chix 17| No. 13 
Triangular on (a, b If x and v are extreme valu 
of the sample, then 
K } X \ ; 2 b 
. i—-7r, 
Vai x+-\ - i O(n 
‘ l6n 
\ 7 
Ei (x—v I ] b—a o(n 
2n 
} wr 
Var xX \ b 2 UT 
Hr) 
See aiso yilocle 
5.11. Rectangular (a—h,a+h 
" 
Dix >} »a—h<x<a+h:;: D(x 0 elsewhere 
2h 
m=a, v=h?/3 1|2 13 
Skewness=0, kurtosis 605 + Chix Si] 
136, 1-144 
Chix sin ht/ht)e*" 1\259 
Special case of Ty pe I] 21142 


rectangular over a,b 


12 


bibliography for 


17) Nos 
15\93 


Il, 


, a 


}+x/a 
If x and Vv are the k™ values from the top and bottom 


of sample 


E(x i h 2h E{4(x-+-y i 
n | . 
2k 
E(x—y)=(1 ) (2h 1\372 
n+] 
. k 1 k | 
Var (x 2h 
rn | n-+-2 
tkh 
Var [2(x+-1 
2(n+1)(n+2 
| k+1 
«KUT ai 
Var (x—y 2h 
I} l I 2 
h h 
~ D(x N (0 ); ~Di Laplace (0, ); 
ti I 
where e=} (max+min), n var (« 6 var (Xx 
) ‘ 
D(z d)26:115 
Moments of max and min VR }} 
D(X) given incorrectly ni10-3:91 
Dix No. 5.16 
) 2low TIx lt] 
FD ia Rectangular(max hh, min h C)IO402 
D(q"™ ranking iten Type I 123390 
Testing against Sill ple unimodal distribution 
420) 11] 
Grouping eorrections C1440 
Transformation to Cauchy Lha\10 
(' R theorem mav not hold j LSS 
\MILE(a hoa Hy max, min Ol L5¢ SL IS 
Best linear estimate ot m and a] diZ IS 
Estimatior d\1? 
Location and sealing closest estimate / 22 
Minimax estimate of ‘‘a d\22 
UMVUE 14\142 
Saves’ theorem 16—-1:110 
Variance of estimates of i gst [0 


Testing 2 €|25 2157 
Critical regions 


Vath. 3:3, 


5.12. Rectangular (0, a 
D(x l/a, Chix e*'t_ 1) /ait 
MGF (x sinh Jat bat L0\38 
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Ct, ha 2r+1 \1O\14 
Cumulants [10\41. 
Var(mean deviation) ™~a?/45n 1217 
Mean difference: lc|28:432. 
Ch(Sx,), C(Sx (9278, 
DiGM a log a/x 12 4e 
( a} 
D (range) =n(n—1)a~"x"-?(a—x 
[4)92, 123, [c|20A:210, 
D( SS log x, )=Type II 16) 
D(x, s) for n=2,3 d\3:128 
LD quotient of ranges Q\46 502 
D (imax, /max, No. 8.63 O71 13 
FD (a x x" C\30:408, |d\9227 
. | als 
\IGE ( log log ); \IGE aresinh yX 
x 
Completeness E1074 
Estimation by order statistics d\2¢ 


(Juasi-range A)28:179 
Best linear estimate d\14:88 
. | 
[ MVI Ie a ( | ) TaN ite 
n 
Sufficient statisties el17:214 
Confidence intervals for a 
Example d\11:209 
Estimation of dispersion: Cf 
UMP test of a l 1s 
See also A\2:48 d\2 6 
5.13. Rectangular 
Dix l QOsxsl 
) low x e~*. Ch log x 
) S log x Type IIT (lk 
| ; , 
1)( =x [x x 
} | 


Irwin-Hall” distribution 
1:73, (2|240, 244, |c|19:234 
VR 7:311, [e}l: 41 
Obtained by finite differences 
10:52, 255. 324 


VR 6 VS 


Shi ppard’s corrections 


M240, (¢ 


Convolutions 


\ 


( 


y 1bN17 
ad ] i 
1) 

‘ 
if 
1: 
) 
X 
fils : 
09. 


VR 12 


O°27¢6 
/ 
; 
Ae 
SA 
‘ 
if 


s sufficient 


[14\142. 


Also for 


I \2 
VR 


No 


v2 


) 


t. 


14: 


Ww 


iD) 
$a 


4) 








D(q™ value from top of sample)=Type I, 


m=1—q/(n+1), v=q(n—q+1)/[(n+1)? (n+2)}: 
[2}218. 

Var (¢)=1/4(n+2): [2|230. 

~D(z): [d|25:636, MR 9:360. 

Mellin transform: [d|19:373. 

Order statistics: [c]24:260, [i|33:214. 


[c]25:417. 


(1-+-t) 


9(eitt 


D(range) = Type I: 


MGF (log log 1/x) 
MGF(aresin x) 


Ratio of two ranges: 


1)/(t®+-4) 
(d)21:112, [x7 :179. 
Z 13:30, [c|20:217. 


[1]485. 


[v]8:69 


Moments of the range: 
C.-R. theorem may not hold: 
[c]41:230 
[d| 13:44. 


(t 32172 


Censored sample: 
Stratified sample: 


Significance levels for X: 


D(GM)=No. 8.12: [wl :73. 
Estimation of center, D(X), D(é [¢|33:126 
D(max,;max,) = No. 8.64: (g|50:1142. 
Two rectangulars added: [n|8-3:74. 
Hypothesis testing: [e}32:321 

) 


See Also: 
418, [y|24:2 


[d|23: 43, le 


, MR 7:310, 


Z 11:218, MR 16:602 


5.14. Correlation 
x’)} n ve] ) \ 4 n-|l gn 
D(x) 
rin >)! 
(2xp)'_., /n+1—1 
; & 7 r( ) [1]398. 
i=0 1! \ 2 
yp?) §O- 1] — x2) i -4 q®-3 cos”! ( x) 
D(x) . r 
a(n—2)! d(px)"~? 1 — p*x? 


[2|342, [10]200 


Special cases n=2,3,4, moments: 


n=4: [u}26:5356. 
C(x): [c]25:71. 
~D(x): [n]1—4:1, [e]38:236 
Ifp=0, D(x)=No. 5.5: (6|314, |4]120, [g]26:129. 


If p=0, D( : yn 2) Student(n—2): 
yl e 


[213 43. 


Transform x=tanh z, p=tanh ¢: [10|200, {c)21: 


358. 
If p=0, D(x*)=Beta [4, 4(n—2)]: (10160, 192. 
Bayes’ distribution of p is No. 5.5(0,1): [3}91, 
(c]41:278. 
Moments: [n]5:3. 


5:2083. lm 16:120, [d|22: 
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Papers dealing with this distribution generally 
(c]10:507, [e]11:328 
Interval estimation: le|7:415 


Confidence limits for p [¢|29:157 


Stratified sampling: [i|36(Suppl.):8? 


See also: [6)15:193, [c]21:164, [c]24:383; [o|3:1, 
Z21:41 
5.15. Multiple Correlation 
y (I—R2)* 
D(x a xk l—x k 
Bis n k) (k ] | 
/n—1l n—1 k—1 
where ¥ K( > > ; R X ) and 
‘ al ly b(b+1 
F(a, b, e, x) Pe a a 
{ Ci | 2 
[14]65, [2|384, [d\3:1.96, [2|127 
If R 0), D xX Beta k, hi k 00 ; Si 
Testing [4|2.9:25 
Another form: 121387. [3|338 
Limiting form when n 1387 
When R=0, D(x) =Snedecor: 10\257 ,2 
Mean, variance: ¢|22:353 
Moments MR 14:189 
More generally: [d\11:6 
See also: le]9:352, |e|10:257 aj92:445, |n|1 
4°04 li 24:199, lz 1:67 137, id L:88 Pro« Roy Soc 
A) 121:654, [u|46:521. 
5.16. Rectangular Mean 
l <), 
D(x >, (—1 x—j 1244 
n—1)!f= 
Generalization: [7}.3:330 
Obtained by finite differences p\l0:52, 255, 324 
Called Irwin-Hall, cf. No. 5:13 d\13:43 
Compare No. 8.47 and No. 8.70: 
See also: London P.O. Res. Re P 13443, A chiv 
der Math. 3:3, Proc. Intl Cong Math 1924) 2:795 
6. Distributions on (0, ~) 
6.1. Type VI 
D(x C(x—a)?"'(x—b)? x >b <b, q>0 
p+q<l 17 |249 
Ifb—0, C : 21.40 
B ] p q, p 
Roots of quadratic In Pearson equation real and 


[11|44 


same sign 


Truncation (39:63, [|40:18. 


[11|83. 


Various constants and an example: 


See also [\d|7:23, [c]23:143, [c]25:379, [17|No 

26, |v|4:167 
6.2. Snedecor (p, q 
y/q Px 1 yx /¢ }(p+a P 
Dix pid x >O: |5|100 
B xp, x] 

“F” distribution 

Mm, V, K3,K4,81,Bs p\IGA745 
Derivation, properties, examples [15|374, 


1S|1 


Area unity if p,q both even: 


163 
MR 12:509 


Obtained as distribution of ratio of two Chi-square 


6110.45, |4)113 


variables 


Ses 
eS where n=}p—1, 
x+k : 
m= (p+q)—1, k=q/p, and e k™ 
| p|10:62. 
D \ xX ) 100 
yx tas 
D(—P*_)=Bip—1, q—1 [4115 
q px 
Therefore called “inverted Beta”’ le]33:73. 
Various properties [d|12:446. 
q 
m LO\TIS 
q 2 
risp+r)l !—I ; 
a; q/p vee 
Pispyrisq 
pq—2q . 
Mods LOWS, 
pq+ =p 
Approximated by normal distribution [d]13:233. 
If x, \ each Snedecor nl n then D(,x/v) 
Snedecor (2n-2, 2n-2 
Testing d\13:37 1 
Used to test multiple correlation coefficient 
10\257 
el | iso LTS d 6:204, ld IS:S2, le eis D0, 
>] 4 | did it J Sor Statist Paris 962262 
6.3. Beta of Second Kind (p, q 
\ p 
Dix m 
B | ( | X q | 
pip—+—aq l ) | ) l 
. | , mode=?—, HM=! 
q—1)*(q—2 q+! q 
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= : 
if r<q, a, pip tn. 
q—1)(q—2 (q—I!I 
D(1/x)—No. 6.3(q,p), D{ =) Beta: 
[1)242, [10]156, 158, 163 (d|25:402. 
p=q, xel: [d|22:418 
Called Fisher’s F: pws 
C(x): | p|7:102 
6.4. Hotelling 
») x" 
D(x) = 41238. 
3 Bi }(p q), +q|(p—1)* x ap 4]23 
p—l 
For q=1, this is positive half of Student distribution, 
hence called generalized Student 
(1)409, |d\2:376. 
D(x? ([c]32:70. 
Mellin transform: \d\19:373 


Percentage points, relation with Chi-square: 
\d\27 :1091 
See also: [10|207, |c|25:399, 30:66, |d\9:235, 
[c]24:480, |c]4:174, [c]24:487, |t|7:82 
6.5. Pareto 
D(x) =p/q(q/x)! (8120. 
p i/p ) ) ) 
ay q, €=2°""q 1\248, |2|142 
p—l 
More generally d 7 26 
Testing, location and dispersion (7 :115 
As Type XI c|39:178, |17|No. 25 
Ranking: (c]24:234, 241, 275 
Double Pareto Ie ndall and Buckland, A Dic- 


tionary of Stat. Terms 


Ne) also 19:17 4, lA’ L499 \¢q LS 37. [7 8:76, 
lf 77 ly 13:30, VRI3: 2 4253 :63 ( R Acad. 
Sc) Paris 23321 L271, {/ P5AY1. / gf. 47 

6.6. Kendall 
ia x—I 
Dix O<— 1 x 0<a<l 
af T(x | 
I ra 
m >\ , 
l ( a 
h / j 211 c} Borel Tanner 
6.7. Inverse Gaussian 
D(x) =exp A(X—)*/2y?x|[A/2ax introduction, 
moments, estimation d\28 :362, 696. 








7. Distributions on (— @, ~) 
7.1. Type VII 
(1 j x? 2) m 
D(x) pier sy m>}#: [2]142 


aB(4, m—}4) 


Estimation : [c]36:412, [t]4:35. 


See also: [c]15:401, [c|36:412, 167 


7.2. Student (» 


l'[}(v+-1)] | y2/,)-$v+1). (5]97 


D(x) (2117, {6] 


I'(4v) (xv)! 
10.6, MR8:161, [14]47, [4]110, [10]186, 
[18]1 162, [17|No. 29. 


[1]3:355, 


Introduction, properties, examples: [15}388. 


v=v/(y-2). [1]239. 
Type VII with m=4 (v+1), a®’=». 
*t”’ distribution. 
1-3-...-(2k 
4)-...+(y 


1)v* 

-2k)’ 

Original paper in which this distribution was dis- 
covered: lc|6:1. 

Ch(x), refs for »=3: [17|No. 28. 

(b]18:212. 


Distribution of the ratio of a Chi-square variable to a 


Ox [1]239, [10)208 


(v—2)(y 


Ch(x): 


normal variable: (1]387, [4]110, [10|187 [n}s: 
102. 
In bivariate normal samples 
o,(n—1)! ; 
D( — (b 8) ) Student (n—1): [1]29.8. 
o,(1 p*)) 
s,(n—2)! , 
D( <q (b ) ) Student (n—2): [1]29.8 
s,.(1— y*)3 
if p=0 
*\} v . ‘ “ 
D| (n—2)! — Student (n—2): (129.7. 
(1 vy)? 
C(x): (/|3:358, [c]25:389, [n|5:109, |[c|37:168. 
D(x?) =Snedecor: (6|217, [4]114. 
Transform to Type II: c|28 :308 
D(X): [n|JS—4:92. 
As vo, Student—»N(0,1): [1]252, [3|101, [a 


113:228, |d\27:783, Proc. A.M.S. 6th Symposium 
in Appl. Math. p. 251. 


Approximations: |d|7:210, |d\9:87, 
[c|34:176. 


variables: 


[d|17:216 
D(log x): 
[c]22:405, [e]23:1. 


(10256. 


Two Student 
Used to test partial correlation: 


~significance levels: |d\14:60. 
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Generalizations: \d|19:406, |d|25:162, [1|34:58 


See also: ld 10:265, |d|18:89, [e|11:37, le\12:89, 


[7 S:6 32. le 332362, In 5290. l¢ os: 271 300. ( 24:56 2 it) 
19133:138., \¢ L4:264, By it. Assoc Math. Table N rd 
fed.) V.1 Pp rrriii, J. Soc. Stat. Paris 92:262, MR18 
834, 24:67, |uj21:482,655, |w)129, [p|8:42 
7.3. Normal Regression Slope 
o7o3(1 »7) 2™—-DP(An 
D(x) : : : 
yrl'|s(n—1))o% J 2 pa ,o.X+ 07x 
] {j7 
] a — 
\ a l p 2\9009, ell 52 
n Od 
For p=0, can use Student distribution to test x 
| Os 
LOWLY. 
Stratified sampling (136296 
7.4. Cauchy (p, q 
] ) ahs 
Dix | 1\246 135, |18)1-—3¢ 
nr p+ (X—q 
Ch(x)=exp [qit—ptt [5|60 
q is the mode and median; there is no mean, nor any 
moment: [5|58 
Quartiles are qtp: [5 |67 
X is not a consistent estimate of q D105 


X is a “density unbiased” estimate of q: [d|25 :400 
There are no sufficient estimators: [3]48 


C.—R. (q), C R.(p), C. 


Distribution of t and F statistics 


R. (p/q): [e|8 2205 
MR13:665 
Mean and variance of 4 (x+y), where x and y are 
respectively the k™ values from the top and bot 


tom of the sample: [1|373 
1 


}(X-+¥) Is not a consistent estimate of q 
See also: ld 17 :2, (d 21:133. 
7.5. Cauchy (1, q 
1) (x) : : 6\117 
nr 1+(x—q) 

q incorrectly asserted to be the mean: pire 
E(x), E(x’), Dix+y [14|43 
C.—R. (q)=2/n [1]4.90, [3]24, [p]7:1469 
Var (£)~7°/4n d|6 
D(x D(x), hence & not consistent 312, (114.90, 

(u|22:702 
D(é): [3)46 
MLE # minimax 16\64 
MLE is solution of Z/2(x q l xi—-q () 


(3124, [p]7:169 


15:194 


Gaging: l¢ 


There is no sufficient estimator: NO16 3127 
p 77162 
There is no UMVUE 51 
Information and estimation e|8:315 
Loss of information 332 
Testing m=m AMI:S83, |d\13:65 
Cauchys added VWR17:863 
See also [b)9:61, [(]20:61, [g]51:641, [d|28:832, 
wits 
7.6. Cauchy (p, 0 
Dix p/m){1/(p Xx Chix exp—pt N27 5 
Reproductive property N27¢ 
Information and estimation le|8:316 
Completeness e\lO:514 
D(x D(x nilO—osl 
Truneated p.p LO\LS 
7.7. Cauchy (1, 0 
Dix l/r) [1/(1+x 5 |167 
Chix exp t ANOS 8 167, 1\246 N243 
17| No. 27, {s]213:718 ISS 
Sample median 1\26 600 
D(x) =Dix 21233247, |w|1:73 
From an example S139, [91242 
Moments SII 
As distribution of ratio of two normal variables 
TO\LS9 
C(x , l/r) tan7'x 
Censored sample ¢)41:250 
Wrapped-up Cauchy d\2 ) 


See also d APS. j 7:371. S DD Poisson 
Is Pi Su la / aret lite de ° re sultats moyens des 
obser ations.”’ Connaissance des Ti ms ou lh S Voure- 
ments ¢ élestes ad l’usaae aes astronomes et des Vavi- 
gateurs pow lan 1s »? Le Bureau de ~ Longitude x. 
Paris d)22:418, MR3:232, MR9:235, MR5:124 
7.8. LaPlace (m, v 
Dix 2\ exp| x—m|/3(2v)}] 1\247. |5|35 
67. |\18|1-—136 
Ch(x exp (mit 1+ vt 62 
Mean and variance of average of greatest and least 


sample values 1|375 


VRS 
and (1 


A priori distributions of m,\ 


“Best eg 


estimates of m and 


rare £ n)X\x iS 
but £ not sufficient 4|147-8 
Distribution of smallest sample value g\43:408 


D £ d ft) 
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Quasi-range A\28:179 


“Double exponential” distribution Proc. Roy 
Soc Lond Series A 154:124 72164 

Convolution, estimation, generalization Sobre 
la Primera Le V de Errores de Laplace, } A Sales 
Vallés. Thesis. Barce lona 1 14,7 

Laplace 0, v), also called ‘Poisson’s first law of 
error,” D(j&}), DOx|—\y Dilog x\—log \y 
D(>x?), D(GM), D(HM d\6:102 

Laplace if D(x)D(v)=¢(\x|+ \y WR10:125 
See also nilO—3:S0 

Laplace m,2), MILE 3 Ai) 

Laplace (0,1), D(x he 81120, |\9\279, |1|100 

Sample median: A\ 262599 

C1, 2k ld 5 

D range quality control PIO 284 
See also (d|22:425 17 Vo aS, Laplace 

1774 Me MW Ere Sul la Probabi lj te des COURSES pal les 


ChEHE METIS 


7.9. Fisher 


P, q 
») I 4, Ps 
2p iPaq is 
DD x : I 1 ] 2 4. ? 249, 
BCP, 34)(4+ pe’ 14\48  MRS:161. 
z distribution: (3:355 
(C(x) as a power series (/|.3:360 
Cumulants: MR9I:48, 735 
Moments, cumulants: [y]5:317 
Transform to Snedecor by ve": 
x i log GX pxs (c]23:147 
Various properties d\12:429, |c\384:173 
' q\'"' 
Ch(x rls ( iwrs(p+it) /rapyrade 
& 3(q pq 
1116. 
~C(x d\28:504 
Obtained from two Chi-square variables [2|249, 
d\7 742, 6173 
Approximate significance levels, transformation to 
normal [d/)11:93 
Normal limit Am. Math. Mor thly, 50:100 
Generalization: ((|34:58 
Noneentral: l0|7 2:57 
See also l¢ 2-423, lh ] i]. a IL: 254 / i¢ 21: 3 0, 
C)d4 352. 359. \¢ A | Am. Math Monthl y ¢ 5O: 100, 
IS2 Pros Intl (ong Vath L920\805. Current 
Scrence 2 p 1a]. lo 5:30. p\6:1ss 
7.10. Type IV 
D(x) =C(1+-x*/a*)-™ exp (—p tan! x/a [11)69. 
D(x h Type LV I3\48, 








Various constants with an example: [11)69. 
Roots of quadratic in Pearson equation complex: 
(11}44. 

a,=|a/(2m—r—1)][(r—1)aa,_.—pa,_;]: [2|86, 140, 

144. 
a=): (17|No. 52. 

See also: [d]7:21, [c]1:39, [e]3:312, [e]6:435 

[e]7:74, [c]26:386, 

8. Miscellaneous Univariate 
8.1. Pearson 

1 dD(x) fig+ a, X 
D(x) dx bo+ b,x+ box? 

types listed with associated parameters:  [d]2:3.94, 

[n]11—4:77. 

Differential equation for Ch(x): MRS :3.93 
Ch(x): MR10:705. 
Original paper: Phil. Trans., 1895. 


First seven types treated in [11]. 
[d\7 316. 
61. 


[d|22:256, |c|40:40 


New classification: 
D(x): — [d]18:111, [n]8-4 


Truncation, estimation: 


Bivariate: MR9:363 452. 
Flexes equidistant from mode, ete.: (d\6:1. 
Bivariate generalization: [v]3:273. 


Orthogonal Polynomials: Ann. Soc. Cien. Argen- 
tina 156:3. 

Generalizations: [d|5:124, [c|26:129, (cf No. 8.52), 

MR10:386, MR17:1095, |d\21:289, J. Gakugei 

Tokushima Univ. Math. 5:29, |n\12-2:95. 


Log Pearson distributions: Intl. Cong. Math. 
(1950)1:4580. 
Romanovsky’s generalization: [e|17:106, [e\18: 
221. 
See also: ld S:18, ld S:206, ld 20:461, l¢ 6:415, 
[e]7:127, [el16:106, [e|16:198, [e|18:264, [ce|20:389, 


[g126(P):288, |a]85:488, |c\35:113, [c]32:81, |c\36:151. 
[e]38:4, [i]25:141, MRI17:169,272, Z9:314, Z6:268. 
Z19:73, MR 4:756 977, Intl. Cong. Math. (1950 
1:585. 
8.2. Bessel Function 
Mahalanobis’ Distribution (“D” distribution l¢ 
2:143,385, \e\3:105, [e\4:19,373,535, [e]8:167, [kIs 
379, [14)246, MR4:23 
Wilk’s distribution of dispersion determinant, ete 
le 13:26. 
D(x): MR5:42. 
Distribution of vector correlation: [c]28 7353 
Distribution of difference of x? variables (d\7:51 
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D(x), Ch(x MRI4: 


nn’ 
‘@ 


Moments in a special case: 
bd. 
Prese nted 


D(ratio Type III variables Studies 


to Richard von Mises. Academic Press, 1954. 
Pp. 301. 
Bivariate Gamma distribution: [e]5:140. 
Distribution of the range: (d|18:384, [d]21:133 
Marginal total of Elfving’s distribution: [c]36:142 


+22 


See also: 


lel6:175, | 2235, [a)97:125, la|98:89, 
(b|16:96, [c]21:168, [c|24:441, [c]24:485,492, [d|18:392, 
[c]24:39, [c]24:293, [17|Nos. 10,60,61,62,63,64, with 
Chir). re fs. lc 21 :164, le 24:293, le 24:39, Several 


Par 18 6: Phy 


forms, J. Soc. Stat 
607, MRI16:152 


2, [uJ28:458, MR11: 


Variance Ratio 


8.3 


2(1—p*)38- ’x®-*(1 tp*x*(1 
B[}(n—1), 4(n—1)](1+x?)" 


vad 
ZIS65, \¢ 


D(x) 


2:65 


‘““Bose’’. 
0, DI 


For p 2n—3)x Hotelling (2n n—!] 


See also: | p|JG:183, 


|p 7:98, [c)380:190, [e\3129 


8.4. Kullbach 


nx"? 


D(x) - 
rin) l 


Distribution of GM from Type TIT (1,p):  [2|25 


8.5. Noncentral Student 


le ‘173.307 


D(x 31:362, [18|1—-162, [v4 
Multivariate (v]4:331 


Application: [c]43:219 


(\36(Suppl.):21, |r|1:28, J. Soe. Stat 


VRI D3 6 


also 


96 2262 


Spe 
Paris 


T 


8.6. Continuous L: 


x1cn 


D(x f(p)()p*(1—p)" “dp; parameters; 


Beta, D(x) is hyvpergeometri NST31, [34219 


8.7. Noncentral Sa-d 


DD xX ( J6: 178 ] 160, ‘ 15: 32] 


AHA) 


see also ( 38:112. N36 Suppl ,. 4 irl 


Fisher's Logarithmic Series 


| x 
‘ 


Dix »X 
log 
P08 Tk 


negative binomial (1, 


Cf m). No. 3.5 0137-358 











8.9. Rank Variate 


} 


‘J 


N n! x 
Dix : , | exp n—q-+ | 
a (q | n—q): 0 2 
| e-~ee 0124:231, 239, \¢ : 
In special case called Yule’s distribution, MLE 
l¢ 42:23 4265, \¢ A3:248. 


A distribution of Type ITT median 


8.10. Generalized Pareto 


p 7 :143 


largest 


function 


similar 


Z4 


veneral distribu 


/ 


D(x)=ax-"{1/(e"*—1)]: [62184 
b=1, n=5, a=15/r', Planck’s radiation 
17|No. 56 
8.11. Ghosh 
‘) 
Dix x Ml where [k] is the 
ilk 
teger * k 
Furnishes counterexample to theorem on 
regions EIS 2390 
8.12. Rectangular Geometric Mean 
nx - 
LD) log a/x 2246, (d\5:276, 
“Ein 
8.13. Cauchy Median 
2m+1)! (a | 
Dix ( tan” '(x—k ) 
mi!)*7 } | x 
PA 
8.14. Spearman’s Rank Correlation 
Dix C130:256 ( SA2183, Clos :1lol 
~Di(x Type Il, moments C\40:409 
8.15. Circular Normal Correlation 
D(x n(n-1) « l-e— x ( 
Ls if 
8.16. Koopman 
Dix \) k) Rix exp k Hix most 
tion admitting a sufficient estimate of k 
d ,) {} ( fy 7 ] Trai \ An Vath Sor 
) li2 
8.17. Von Mises 
L(x C exp k cos (x—a Physikal 
PO:490 
Circular Normal GQi4S:1ol, |q\4l 26 
a) d\26:233 C43 :344 
8.18 
Family of distributions having all moments equal 
| 
J] {j? 
Dix ( lI—psm\ O<x (} 
tk 3)! 
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8.19. 


Family of distributions having all moments equal 


D (x) =e7! w) x7! * 1] » sin (2m log x 
| | 
Oy exp Lk k+2 id 11:402 
8.20. 


distribution”, involving a hypergeo- 


“Non-null t’ 


metric function: [14|48 
8.21 
sech*~*x , ; 
D(x : » the distribution of tanh™'r in 
B(1,k—2) 
samples from a bivariate normal distribution with 
zero means and zero correlation 1b)15:213, 
|b GA] 
Ch(x), also special cases and refs [17|Nos. 53-8. 
D(2x,) for k=3: ZI:219 
8.22 
I x—m 
D(x)=< sech? ( ), connected with lognormal: 
ail a 
é 12 ] ,? 
Called ‘logistic’, Ch(x rx sechax when m=0, a=2:; 
V7 21638 
Called ‘hyperbolic law of errors’, properties, ex- 
amples 8|2:55 


UR 13:665. 


Distribution of t and F statistics 


8.23 
Distribution of the correlation ratio, involving 
Series la 17-2121, ( 24:441 
8.24 
Various distributions of the form exp quartic 
polynomial [d4:1, [di4:79, [d\19:589 


Gri\ ge an example where ho minimum variance esti- 


nator exists €|12:43 
MILE CWTTSS, |aj\IS:114 
8.25 
k(1+-x 2|52 
Ch X 2104 
8.26 
D(x b sin 2(a bx b\9:61 
8.27 
Normal multiplied by an eighth degree polynomia 
a\lOog-s6] 
8.28 
1) X th th X }¢ nix 1D) x nilO—S 2:90 
' ; | " ' 
For h=1, Ch(x)=(, ) 10:75, (17\No. 39 





8.29. Determination of constants: C.R. Acad. Sei 
Paris 222:34. 


Miscellaneous distributions given in terms of C(x): Esti MLE 124-989 
(d\13:217, Math. Tables and Other Aids to Computa- SUMAWON, SLE 12) 4-208. 
tion 5:109, |g\50:209. Bias: [d]27:758. 
8.30 Ch(x) =(1—ixa™') exp(ixm) Cumulants (see refer- 
— ence): (78]1-144. 
D(x) = $k(1+- |x|)“*"!: [1]225 No. 2. For m=0, a=1, Ch(x r(i—ix), References 
. (171No. £3 | JS. 
k negative: (c]33:126. (17|No. 43, [w]24:180. 
k=1- p: [c]36:93, [17] No. 17, generalization No Connection with No 2.3: lv]4:8, |glo0:518 
18 : [9]42:408 
8.31. Special cases Ir|1:4 
D(x) =k exp (—ax”)/1+q sin (bx?)]: [2|106 ~D(x [uJ24:180, |g]43:408, Studies presented to 
richard von A1ises. Academic Press, 1954 p. 346. 
ers ) Richard von M icademic Press, 1954 } 
8.32. With slight modification ~ D(x) is distribution of log 
ne — ; ; : survival time: J. iene 42:328 
Various distributions formed from rational functions ., aaa Hyg sete 
of x, rational functions multiplied by e~'”, e~* and see also: US Dept. Agric. ARS 41:13, Ani 
exp (—tan™'x): [d}1:137. Inst. Henri Poincaré 4:115, 5:115, J. de Physique 
Serie 7? Vol o, OSs, Ss, I]. Bull Am Mi teol Soc. 
8.33. 23:95, CR. Acad. Sci. Paris 246 :49, 237: 512, Nature 
175:270, Cong. Intl. Math. 1936, 2:200, |a\99:732, 
D(x) = (e?+ |x ' {log (e?+-|x|)]~?, having a pathologi- | [w|8:97 NBS Appl. Math. Ser. No. 33 
cally long tail: [d\17:11. 
8.41. 
8.34. Weibull 
2 x? Ch(x) =e" (1 x), [n]L0:75 
D(x)=ab x"! exp(—ax’), x>0, a>d0, b>1: | Dix ry 
[9143:408, [17|No. 44. 
Moments of order statistics: (d\26:330. [17]No.30. 
8.42. 
8.35. 
D XT YV); where x and \ obey Various trivial dis- 
D(x) =k sin™x cos"x: [c]30:182. tributions [|d]5:16. 
n=0, value of k, Ch(cos x): Phil. Magq. Ser. 7, 
39:70. 8.43. 
6.36. D(x) =Cx7'(1+-p/x) [d]6:106. 
D(x) =2hx7!(1+h?x?)~?: [n|10—3:77 
8.44. 
a+ 1 — = 
6.37. D(x) ="— (1—|x|*), —1<x<1: (17]No.15 
ail 
D(x) =2hr'(e™+e-"*)-!, D(x) called ‘‘Perks’”’: 8.45 
[n|10—3:90, (d|26:153, [v]7:159, J. Inst. Actuar. 40. 
63:12. 2 ! n e)x?J 
8.38. D(x) =-————_—-. e~*x® > ( ) ——— » O<X , 
(1 c ) ) / T(bj+a 
D(x) = (2%) ~4(x/k) [exp(— $(x—k)?) »” 
x (xy +-k)?2)]: 21387 
exp( 2\X k) J ¥)000~ X 0, c>0, b -O a oO. n 1.2 
8.39. Chix), References 17|No. 3¢ 
Four distributions formed by multiplying the normal 8.46. 
distribution by a polynomial, used to illustrate : 
kurtosis: (g|40:259. D(x) =(x—k)x" e™ g\42:57 2. 
8.40. Extreme Value 8.47. Stevens-Fisher 
fn n [.19.@172 , 
D(x)=a exp|—a(x—m)]| exp[—exp—a(x—m)]: D(x =(. ) 1)4(1—)x [A]9:315, [k]10:14. 


[d]17:299. 
Gumbel or Fisher-Tippett distribution. (‘ompare No. 5.16 and No. 8.70. 
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8.49 
Dix (* exp i(b—x d\25 
8.50 
My 
Dix = ( Xx 
Or (ia 


26(Suppl.-H )227,{17|No.40 
MLE u)45:542 
8.51 


Distribution of non-normal correlation 


D(x c[P = | \P<x<yp 


Value of .;. relerences 17|No 


Hansmann’s distributions, obtained from 


ized Pearson differential equation 


8.53. 
D(x k/x) exp[—ax— (b/x 17|No.48 
Called ‘Type Harmonique”’ R. Acad Sei 
Paris 213:634 
8.54 
D(x uv yre"’ ( x j23:101 
\ 
Ch(x exp[2a(yb—yb—it)] 17|.No.50 
8.55 
a ’ 
D(x : enn O<—X 1 >0, p>0 [17] 
I'(p 
No Jl 
8.56 
D(x oe ae". i<is loa 17 
Vi 
8.57 
D(x Lb x* lo x] b x vt 1. b l 
and 
ly , 
L(u >u>o 17|No.58 
log \ 
8.58 
D(x i+1)*x*logx, O<x<l, a 17] 
Vo ; 
8.59 


\ veneralization of the livpergeometric distribution 
based on the Whittaker function 
‘ yy 1 


x e F.(m+4—k, 2m+x; 


D(x), Ch(x), references 17|.No.t ). 


vil 


i creneral- 


c]26:129 


8.60 


Dix)=C eW* /(b?+-x*), moments, Ch(x), limiting 


cases Cauchy, ormal r2:293, r3:139. 
8.61 
D(x l/r l COS X)/X 22020 
8.62 
D(x) where A+ Bf(x) is (a) normal, or (b Laplace 
and f(x) is (i) log x, (11) log [x/(1 —x)], (iii) aresinh 
x [c]36:149, |r]5:283 
8.63 


8.64 


Garwood’s distribution of length of gaps in traffic: 
b\7:65, |g|46:117, |c]38:384 


8.65. Matching 


D(x 1/x!)[1 1/1!)+ (1/2! ee 
1)*-*/(n—x)!] aj118:390, MRI18:346. 


Generalization 4\210 
8.66 

Cigarette card distribution a\118:391. 
8.67. 


Cubic polynomial over a finite range, estimation: 
ga0:196, d| IG FOS JG] 


8.68. 
mn 
0< x< 1 
m-+-n 
D(x 
nn 
Dp : l= x 
m+n 
\loments. ete g\b0 1137 
8.69. 
mn(m—n)~'x' |—x’ »> msn 
Dix 
Ln x""' Jog (1/x), m=n 
\loments g\o0:1142 
8.70 
n' n } 
Dix » De l x »Os<xs1 
n—1)! jn ( ) )( ) - 
(‘ompare Nos. 5.16 and 8.47 26:713 








(r—1)!(r—x) _ 
xiy°"* win 


D(x)= 0,1,..., (r—1) 


[6] 19:339. 
8.72. Arfwedson 


Moments, approximations: 


— . x ; 
D(x)=30 (x—j)"(—1)'(. 34:121. 
2 (x—j ( ) tis 
8.73. Stevens-Craig 
D(x)=C n™’o?, of being Stirling’s number of second 
kind: [k]8:57, [c]40:178. 
Generalization : [w]7:203. 
8.74. 
D(x)=}3+x', —1<x<l: [m]6:120. 
8.75. Ising-Stevens 


m 


as...) 


D(x [k]9:10, [d]11:370, 
—s (]4:171, [w]8:55. 
m 
8.76. 
] 
D(x Ann. of Math. 27:18. 
mym*— x? 
8.77. 
D(x) exp [ (x! b)*] Ann. of Math. 
Sy 2r : 
oa ae: 
8.78. Negative Hypergeometric 
n\ B(p+x,q+n—X) 
V(x : 
D(x) (.) Bip, q (6]10:257, Proe. 


Int'l. Stat. Conf. Rome 1958 paper 71. 


Obtained by assuming binomial probability to obey 


Beta. 


8.79. 
D(x)=k(1+x)~? over various ranges: d|22:425. 
8.80. 
D(x) = (n—1)[{1—(1/x)®~*](1/x*), x>1: [d]22:425 
8.81. 
D(x) log x*/x*(1—x°*), mi x<_¢ {d] 22:425. 
8.82. 
a~—bx)*/2ex 
D(x in 
y2rex® Called “inverse Gaussian 
Nature 154:453, {uj43:41, Virginia J. Sez. 
p24 ‘ 


(new series) 7:160. 


Zeit. f. Ph ys ik 31:258, 
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Written 
C(x m)* 
exp ( 2m*x 
Dix : \d]28:362, 696 
(= ) 
D(x)=k xte-****, called “Halphen,”: Publ. 
Inst. Statist. Univ. Paris 4:38. 
8.83. 
D(x)=a(a,p) lexp(a*x*) — 8]: MR16:381. 
8.84. 
D(x)=a exp [—k? log? p] where 
p=|(X—Xo)(X2—x;)]/[(x2.—x) (x, —x Z10:3138. 
8.85. 
(Ch(x)=1/cosh t, t/sinh t, 1/cosh? t: MR11:443 
8.86. 
D(x) =2A(A+1)/(x+A—1)(x+A)(x+A+1), A 0) 
x=1,2 m=1-+A, v le]18:353 
8.87 
N!n! 
Dx [s]5:161 
IT x,'(i!)* 
8.88 


D(x)=ce™/(1+e)* | 


Hy perbolic Error distribution. 
\ discrete distribution from an urn model 


d| 22: AS) 2 


[w]7 71738. 


8.89 
2 l 
D(x)= > mm AXae n' 0<x< oe, 
ax ; '(n 
] ] v-+-a 
O0<a<1, A>0, m \ 
A(] x \ y 
8.90 
x— 1 
—_ p 
Poy. 
Dix)=k*¥ > : x=0, 1 
x—)j)! 
Called “Pélya- \epphi”’ €|40:206 
8.91. 
| —x—2 
D 7 k* 1-+k)** 2 
. A x | ) . 
8.92. 
(/ x 
1 (5) (er (N, A+p)x)/P (N 0 
D(x , 
IN 
| (4 ) (us _ 
L he 


D(vy—z) where D(y)=Type IIT (A, N), D(z)= Type 
ITT (u, 1 
x /fi+N—1 
Or, -*k! Ss . ) (—u/r 
tN 
«= (k—1 ‘4 (—1)*/y 
8.93. Miller 
m+k+1)! (m+x—1)! 
D X K ] 0.1 
m m+k+x+ 1)! . 
m-—+-k m+k+1)(k 1)(m ] 
o 
; k k?(k—1 


Obtained by ussumMinge Pascal paramete rto be Beta 


9. Miscellaneous Bivariate 


9.1. Cauchy Bivariate 


Dix, y wr) [1/(1+x?+ y } 
1 
Dix, v —(1-+x?+y 8 :235 
9.2. Student Bivariate 
x 2rxy \ 
Dix, vy Ix |—) ( | ) 
: nil ! 
Tables ¢]22:408, [e)41:154 


If x,v independent 


9.3. Poisson Bivariate 


Discussion A156, {7:41 f C\o. 146, Ps jt h 
Bull Afi 4 a4 Prog kd } Vii th. Soc ITs 4:18 
Special case obtained from binomial |17 Le 
9.4 
Lognormal bivariate ¢]22:130 30) 
g5 
Normal lognormal 4 O) 
9.6. Binomial Bivariate 
a ef ih ¢ Xp Is { bye ce dl 21133 
If 0, D \ k! i¥(] — 
X | ) ( 
x! y! k x y vii | | | 
etc 17: ij? 19:209 
See also 36:74, MR14:995, Z18:154, MR13: 
BAS 


9.7. Gamma Bivariate 


[e]5:1 40, 


9.8. Gamma-Normal 
l—s exp| it ( ] = )| elé:144 


9.9. Hypergeometric Bivariate 


5 oy eK) 
OE i us 


Cc 16:17 ‘AR c|22:140 


MGF (x, y 


Cc P52] $2? 


rious propert 1es: 


17 :104 


Moments Ganiia 5:97. Koninkl 


A) GO:124 


9.10. Negative Binomial Bivariate 


Various properties 117 :100. 
| (x+y ) m ; 
D(x, y yr 1 " 
p+2m)p [(p) P(ix+DP(v+1) EL p+2m 
m : 
correlation reeression ete c\41:79 


MRI? :605., 


Polva Ke r¢ nbe revel 


9.11. Elfving 


1) x \ ) X eXp Xx cost y . connected with ~D 
range c)34:111, | 36:142. 
9.12 
Dix.\ C es l1—x—\ l+x—y)* ell 4: 
Rhodes surface 22:13 C|41:550 
9.13. 
D(x, 1+x/a)' 1+v/b x+v)/ce, Filon- 
Isserlis surface €\15:222, \e|16:180 
9.14. 
D(x. xv)*¥ (x—v)[(1—x) (l—\ C)O1 2226, 
]-] 
9.15. 
l—y =< 
D(x, IN3665 
|—2rxy+)y 
9.16 
n! 
Dix, " </I v/n)*""(1 
, k 1)']7(n—2k)! 
x/n—v/n)"-* 
x >0, y>0, x+y<cn, 2k<in 
As no, x,v-—>independent Type LII(1,k 


j 
Spec inl case, K 0) Gli i147 











9.17. | 9.32. 
Dix,y)=C xt®* (y—x)i*®"” exp [—4(k—1)y], | Defined over (0,0) (0,1 1,0) from urn model 
D(x/v)= Beta: [b]7:213. [7]3:328. 
9.18. 9.33 
Uniform bivariate, triangular bivariate: [c]24: | “Correlation by common factor” surface c\24: 
382, [v]5:322. 288 , 
9.19. 9.34 
Gram-Charlier bivariate: [e]36:177. Johnson’s system; ten surfaces obtained by trans- 
lation [c]36:297. 
9.20. 
The fifteen constant surface, 9.35 
(quartic lal). e~®@: 117 :268. 
its ; polynomial le] 17 Nine surtaces with Pearson or Jesse margina 
9.21. distributions VWR5:126 
Pearson’s Student-like surfaces: [e]15:234, lel18: 9.36 
229, (c]22:137. 
D(x,v)=C|(x—1)"(h—x) "(v—1) "k—v)! Hoel, 
9.22. Intro. to Math. Stat. 180 
D(x.y)=x+y, D(x+y): [S]94. 
9.37 
9.23. 
: Type IIL bivariate, with discussion and calculation 
Normal-negative binomial: [k]13:289 of Di el? :159 
9.24. 9.38 
Edgeworth surface: [c]38:220, \e]17:314. 
[c]3 Ni ’ Dix.\ '(1+kxyv),|ki <1, 1<x,\ | rs 
9.25. 
Rayleigh bivariate Electrical Engineering, No- 10. Miscellaneous Multivariate 
vember 1954, p. 1004. 
9.26. 10.1. Wishart Trivariate 
, i ta Sd , x Z 
Discussion of “possible’’ bivariate distributions, | D(x,y,z ~ : ayy exp =r V.X—2y 
Narumi’s system; [c]15:77, 209, 222 in (n—2)M 2M 
Generalization (c]22:109. V5 ); 
9.27. Von Mises-Fisher Distribution u= pe M=v,v.(1—, 111397 | 3|3d0 pat 
Generalization of No. 8.17 Proc. Roy. Soe Th ull 
Lond. Ser. A 217 :295, |c|43:344 \ 2M 2M 
9.28 Ch(x,v,z ( .) where A ind 
\ 1 Vion 
Beta bivariate [7]2:261 2M 2N 
9.29 nv ul , 
is —j 
il 2NI 2M 
e7! 
D(x,v) =k [v]3:153 \ 
m?* + (ax?+ 2bxy+ey* ' _ wel 
il 1u 
9.30 2NI 2N 
D(x, } k i] aunxX b*y 2abrxy]", n l 0, 1 l Moments and cumulants 133 
[r]3 273. 
As distribution of normal bivariate variance-covat 
9.3. mines clo a] } e|21 ]f 4 ¢|27 AGU) 
D(x,y) =k exp [—Q(x,y)]-[h?’+Q(x,y)]": [r]3:273 See also k}9:243, |u) 44:295, J. Soe. Stat. Par 
» 96 221 7) u\29 PEL. aj92 5S) 


[r]5:323. 
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10.2. Wishart Multivariate 


Dix a” XK @-* exp Day ;X 
where X xX ; \ a ind 


K,.—**- (Pn 1 r ((n—k 
, A 

Chix (\s) 

where A* a Lé;,l and « yt 


[4]226, [3]331, [1] 391-4, [14]66, [1]30:151, [u]29:260, 


27 I 


Reproductive property 4\232 

Various properties ¢€)20:32, |124:185 

Non-central Wishart Proce. Roy. Soe. Lond. Ser 
See also A3:197, |d|15:345, |d|17:409, |d|19 
42. 13:25, [a]97:120, lel24:476, [136:59, [RIE2244 


COR 0) / ln 440, 
c)38:470, i)86:17, [u]35:336, MR10:387, {b]17:79 


E(x np; Var(x np,1l—p N52, [14]35 
Moments Bull. Amer. Math. Soc. 41:85? 
Introductory article wit! applicati ns 14\|3¢ 


PGE 18\1—-146 

(‘hi-square test e113:2, |e]386:118, 
Information and estimation els 
MLE E\ 182139 


Distinguishing between two multinomials, asymp 


totic torm 14 240 

Trivariat 1146 21 :420 

Biv rinte multinomial }2 Le kad 
de Ciencia Lisboa 2 S 1) 2:19 


Gramma multivariat 


10.6. Student Multivariate 


Student for p=1: [e}41:153, MR16:602 
See also [w]9:143 
10.7. Cauchy Multivariate 
D(x,. , x,)=C(a*+-x?) -* where x?=2x?} 
c]41:54, MR16:51 


D(x,, a ° C, ( l+>° ) : [w)8:235 


10.8. Spherical 


where I y DXi p yz €\41:40 


10.9. Poisson Multivariate 
Derivation e\ll 120, d\28 L060 137 ] 
Without correlation, multiple Poisson 


D(x, a 


‘\ kee 
exp k, k, , ; ; 
\ Xy 


(7127, |e]19:210, Z12:113, 410 


10.10. Binomial Multivariate 


MIGE (x - 1+ p,t Vpyytyt N 
e]11:119, Z12:118, 410 
Dix), ete., in special case 178:271 
10.11 
Negative binomial multivariat« 178:274, [2]19: 


l] Kor j Vederl Akad 1)60:121. 


10.12 

Multinomial multivariate ¢) 36:47 
10.13 

Hotelling multivariate hk} 9:258, |x]7:70 
10.14 


\lultivariate distributions obtained from the normal 


multivariate 127 2 ISS AD 
10.15 
Generalization of No. 9.14 k\9:245 
10.16 
Generalization of No. 8.3 helt 7:13¢ 








Generalization of No. 8.6, No. 


Hypergeometric 


10.17. 


9.29 


- [7]3:153. 
10.18. 


multivariate: [e]15:391, [f\13: 


488, MR17:634, MR12:722. 


Gram-Charlier multivariate: 


D(x): 


Tolerance limits: 


10.19. 

MR14:486. 
10.20. Run Length 

[d]11:367, |4]202, 206, |s]5:1 


~~ 
~ 


Beta Multivariate 


(4194. 


10.21. 
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Selected Bibliography of Statistical Literature, 1930 to 
1957: IV. Markov Chains and Stochastic Processes 


Lola S. Deming and D. Gupta* 


September 29, 1960 


rhis e fourth 1 i series Of bibliographies dealing with various specif 
Viarkov el ns 1 stocha c¢ processes have been take from technical journals pul 
1930 in the t countris hat | e been active engaged tatistical . 

Two prominent reviewing journals Zentralblatt fin the Statistical Engineering Laboratory who origi- 
Mathematik (1930 to 1940) and Mathematical Reviews | nally conceived the plans for this project, continues 
1940 to 1957 have provided the soure inaterial with encouragement and assistance 
for this series of bibliographies. The NBS Statistical The references given here contain the following 
Engineering Laboratory maintains a current file of information taken directly from the abstracts 
abstracts cut from these journals and pasted onto Autho The author’s surname, followed by initials 
cards, on thi subject of probability and n ithematical only In the case of multiple iuthorships the 
Statistics Dr h ecard is coded by subject after the journal reference appears with each author’s name, 
classification svstem of \athematical Reviews, and but the title of the paper appears with the senor 
filed by author \uthors’ names and references are tuthor only The svmbol @ preceding thi surname 
transcribed to punched cards, and classified by | denotes multiple authorship 
subject code to facilitate use of the abstract cards Title: Exactly as in the reviewing journal. Titles 
The list ngs therefrom, tor anv chosen subject, when ol se parately bound publication = books reports, 
supplied with titles, provide the manuscript for thes theses, etc.) are in italics, followed by the publisher 
published bibliographies The bibliography pre Reference to literature Some very inconventional 

sented here is on Markov chains and. stochastt abbreviations were necessitated by the limited space 
processes lor transeription to the punched cards The name 

Complet coverage ot anv one of the subjects is of the journal in italics and the number of the volume 

of course, impossible, but the reviewing journals in bold face, are followed by initial page number 
have undoubtedly abstracted the most important Date of publication: The next figure, in parentheses 
publications fron a large variety of technical journals shows the date when the article or book itself 
and publishing houses. The reviewing journals | appeared 

have vw nerally Invored thre p iblic tions i) Statistics \V/ lol Vathemat cal Be ews ma lor Ze tral- 
that deal with theory und 1 ethodology anv sub hlatt f Vathematih ur followed DY the volume 
Classification on applicat is consequently very number and page number of the reviewing journal 
limited in which the abstract appears 

Following the subject classification of the athe 
matical Be ews Al dal Indes all abstracts that 

come under the general topic of Markov chains and | A. Stochastic Processes—-General Theory 
stochastic processes were consolidated, until 1956 

Into one category Since that time, however, the @Aczél, J See JAnoss) Acta Matl Acad. Sei 
inde. has subdivided this reneral group into \ Tlu t/ 2. 209 1950 

Stochasti processes eeneral theory B Niarkov Adhikari, B. P., Quelques proprietes des processus 
processes c* Stationary processes 1) Special stochastiques localement continus en probabilité, 
processes nel owe introduced = thre fiitt itegory CR. Acad. Se Paris 244, L000 (1957 AY 1s, 944 


KE Appli ations to include whateve publications @ Ambrose W.. ser Doob Ani Vath 11. 7a 


on application Ss that the reviewing journals chos to 1940 

abstract Ambros W On measurable stochasti processes, 
The task of subeclassifving the accumulated ab Trans Ame Vath. Soc. 47. 66 (1940 ML 1, 149 

stracts from 1930 to 1956 into these fine groupings (Anderson | \ A note on a maximum-likelihood 

was accomplished at the Cathole Unive rsitv of estimate, / mometrica 15, 241 (1947 M 8, 593 

America under the general supervision of Professor @Anderson, T. W., Estimation of thi parameters of 

Eugene Lukacs. Dr. Churchill Eisenhart, Chief of a single equation in a complete system of stochastic 

- equations Ani Vath Statist 20. tf} 1949 

*D ( 1 ‘ ‘ 
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